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ABSTRACT 


Hydromagnetic  phenomena  are  induced  in  the  ocean  environ- 
ment by  motions  of  the  electrically  conducting  seawater  across 
ambient  magnetic  fields.  This  report  surveys  and  estimates  pos- 
sible effects  from  a fundamental  theoretic  standpoint.  The  ocean 
motions  considered  include  those  of  natural  origin  and  those  asso- 
ciated with  moving  ocean  vehicles.  The  ambient  magnetic  field  is 
that  due  to  the  earth,  plus  that  due  to  the  vehicle  if  it  is  mag- 
netizable. These  various  motions  and  fields  produce  a broad  multi- 
tude of  ocean  magnetic  (OM)  phenomena.  The  manifold  of  OM  effects 
is  surveyed  theoretically,  and  from  basic  physical  principles  we 
derive  simple  numerical  estimators  for  the  effects.  The  poten- 
tially more  significant  include  several  pure-induction  types  (non- 
propagating) and  several  pseudowave  types  (OM  induction  effects 
propagated  by  mechanical  waves).  Detailed  analyses  of  these  OM 
phenomena  will  be  presented  in  subsequent  reports. 
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1.  INTRODUCTION 


Seawater  has  moderate  electrical  conductivity  so  that  ocean 
motions  across  ambient  magnetic  fields  induce  electrical  currents 
and  associated  magnetic  fields,  electric  fields,  and  charges.  Such 
hydromagnetic  (HM)  effects  have  long  been  known,  are  firmly  baaed 
in  electromagnetic  (EM)  theory,  and  are  readily  measurable  (e.g.. 
Refs.  1 through  4).  Nonetheless,  because  the  oceanic  electrical 
conductivity  is  relatively  small,  the  effects  are  far  weaker  than 
ordinary  laboratory  HM  or  MHD  (magnetohydrodynamic)  phenomena.  To 
emphasize  this,  we  use  the  term  "ocean  magnetics"  (OM)  for  the 
phenomena  to  be  considered  here.  These  OM  phenomena  are  distin- 
guished from  the  oceanic  induction  often  considered  in  geomagnetism 
(cf.  Ref.  5)  that  arises  from  time  variation  of  the  earth's  field 
rather  than  from  seawater  movements. 

The  OM  phenomena  treated  up  to  now  have  been  those  arising 
from  the  interaction  of  natural  motions  of  the  ocean  with  the  geo- 
magnetic field.  For  example.  Weaver  (Ref.  6)  has  discussed  wind 
waves  and  ocean  swell,  Beal  and  Weaver  (Ref.  7)  internal  waves,  and 


Ref.  1.  K.  C.  Maclure,  R.  A.  Hafer,  and  J.  T.  Weaver,  "Mag- 
netic variations  produced  by  ocean  swell,"  Nature,  Vol.  204,  1964, 
p.  1290. 

Ref.  2.  Ye.  M.  Groskaya,  R.  G.  Skrynnikov,  and  G.  V. 

Sokolov,  "Magnetic  field  variations  induced  by  the  motion  of  sea 
waves  in  shallow  water,"  Geomagn.  Aeron.,  Vol.  12,  1972,  p.  131 
(translation) . 

Ref.  3.  N.  A.  Cochrane  and  S.  P.  Srivastava,  "Tidal  influ- 
ence on  electric  and  magnetic  fields  recorded  at  coastal  sites  in 
Nova  Scotia,  Canada,"  J.  Atmos.  Terr.  Phys.,  Vol.  36,  1974,  p.  49. 

Ref.  4.  M.  Klein,  P.  Louvet,  and  P.  Morat,  "Measurement  of 
electromagnetic  effects  generated  by  swell,"  Physics  of  the  Earth 
and  Planetary  Interiors,  Vol.  10,  1975,  p.  49. 

Ref.  5.  E.  C.  Bullard  and  R.  L.  Parker,  "Electromagnetic 
induction  in  the  ocean,"  Chapt.  18,  Vol.  4,  The  Sea,  Wiley-Inter- 
science.  New  York,  NY,  1970. 

Ref.  6.  J.  T.  Weaver,  "Magnetic  variations  associated  with 
ocean  waves  and  swell,"  J.  Geophys.  Res.,  Vol.  70,  1965,  p.  1921. 

Ref.  7.  H.  T.  Beal  and  J.  T.  Weaver,  "Calculations  of  mag- 
netic variations  induced  by  internal  ocean  waves,"  J.  Geophys.  Res., 
Vol.  75,  1970,  p.  6846. 
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recently  Podney  (Ref.  8)  has  generalized  these  analyses,  while 
Larsen  (Ref.  9)  has  analyzed  tidal  wave  effects.  A recent  review 
and  bibliography  of  natural-motion  OM  studies  will  be  found  under 
"hydromagnet ic  noises"  in  Akindonov  et  al.  (Ref.  10). 

Further  OM  phenomena  arise  when  a moving  body  is  introduced 
into  the  ocean  environment.  On  the  one  hand,  t’  * ocean  motions  in- 
duced bj  the  vehicle  interact  with  the  geomagnetic  field  to  produce 
OM  fields,  currents,  and  charges.  On  the  other  hand,  when  the  ve- 
hicle is  significantly  magnetized,  as  the  ferromagnetic  construc- 
tion of  ships  and  boats  generally  permits  via  the  geomagnetic  field, 
the  interaction  of  the  body  magnetic  field  with  ocean  motions  can 
produce  further  OM  effects. 

We  plan  in  this  report  to  survey,  summarize,  and  estimate 
the  broad  spectrum  of  OM  effects  that  can  arise  from  the  interac- 
tions of  either  natural  or  body-induced  ocean  motions  with  either 
geomagnetic  or  body-magnetization  fields.  Our  approach,  through 
theoretical  fundamentals,  recommends  itself  because  the  basic 
physical  principles  of  the  problem  are  well  established  and  mathe- 
matically formulated.  We  begin  by  reviewing  the  basic  governing 
equations  (Sec.  2)  from  which  the  general  character  of  OM  phenomena 
and  a summary  catalog  of  effects  are  outlined  (Sec.  3),  Theoreti- 
cal methods  of  analysis  are  then  presented  (Sec.  4 and  Appendix  A) 
by  means  of  which  we  deduce  simple  numerical  estimates  for  the 
various  OM  fields  induced  by  natural  and  body-related  sea  motions 
(Secs.  5 and  6) , 

In  keeping  with  the  survey  purpose  of  this  report,  the  esti- 
mates of  Sec.  5 and  6 are  designed  for  ready  comprehension  at  the 
expense  of  full  precision,  in  order  to  give  a feeling  for  the  var- 
ious phenomena.  Where  available,  precise  analyses  will  be  cited. 
Finally,  a summary  epilogue  of  the  survey  is  giver,  in  Sec.  7,  along 
with  a glossary  guide  to  the  detailed  estimates,  giving  the  major 
symbols  used. 


Ref.  8.  W.  Podney,  "Electromagnetic  fields  generated  by 
ocean  waves,"  J.  Geophys.  Res.,  Vol.  80,  1975,  p.  2977. 

Ref.  9.  J.  C.  Larsen,  "The  electromagnetic  field  of  long 
and  intermediate  water  waves,"  J.  Marine  Res.,  Vo].  29,  1971,  p.  28. 

Ref.  10.  V.  V.  Akindonov,  V.  I.  Naryshkin,  and  A.  M. 

Ayazantsev,  "Electromagnetic  waves  in  sea  water  (Review),"  Radio-  > 

teknika  i Elektronika,  Vol.  21,  No.  1,  English  translation  in 
Radioeng.  Electron.  Phys.,  Vol.  21,  No.  1,  May  1976. 

»'  t 
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2.  FUNDAMENTALS 


The  EM  field  generated  by  the  ocean  motions  is  governed  by 
Maxwell's  equations.  The  ocean  flowfield  conversely  is  subject  to 
EM  body  forces  through  Navier-Stokes ' equations.  The  combined 
Maxwell-Navier-Stokes  fields  at  sufficiently  slow  rates  of  varia- 
tion (see  Eqs.  3 and  4,  below)  comprise  the  diverse  phenomena  of 
MHD  of  both  cosmic  and  laboratory  interest  (e.g.,  Roberts  in  Ref. 
11).  The  case  of  the  ocean  falls  between  these  two  extremes,  with 
neither  conductivity  nor  flow-scale  nor  ambient  magnetism  particu- 
larly large,  so  that  the  flow- induced  Maxwell  field  is  small  and 
the  EM  force  on  the  flowfield  is  extremely  weak.  Thus  we  have  to 
deal  with  a degenerate  variant  of  MHD,  which  we  call  ocean  magne- 
tics to  emphasize  its  special  characteristics. 


2.1  THE  MAGNETIC  FIELD  IN  THE  OCEAN 

The  primary  relations  for  OM  will  be  the  Maxwell  equations. 


3B 

curl  E + ~ * 0,  div  B = 0 
~ dr  >•- 


(la, lb) 


3D 

curl  H - - j , div  D = p 

CL  /v  ^ C 


(lc.ld) 


where  pe 


is  charge  density,  j the  current  density,  and  E,  B,  D,  H 


the  EM  field.  Equations  1 are  invariant  with  r-ference  frame,  but, 
for  a stationary  observer  of  the  moving  ocean,  its  constitutive  re- 
lations among  B,  H,  D,  E,  and  j must  be  transformed  from  those  of 

/V  «V  +*+  o# 

seawater  at  rest.  Assuming  the  relations  to  be  linear  and  iso- 
tropic with  dielectric  constant  e,  permeability  p,  and  conductivity 
o,  and  denoting  the  ocean  flowfield  as  u,  one  finds  (Ref.  12,  Part 
IV) 


j = o(E  + u x B)  + p u 

M ^ /V  6 ^ 


(2a) 


Ref.  11.  P.  H.  Roberts,  An  Introduction  to  Magnetohydro- 
dynamics , American  Elsevier,  New  York,  NY,  1967. 

Ref.  12.  A.  Sommerfeld,  Electrodynamics , Academic  Press,  New 
York,  NY,  1952. 
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= eE  + (ep  - eftpn)u  * H » 

(2b) 

“ 

* pH  - (ep  - enPn)u  x E , 

XV  U U XV  XV 

(2c) 

* 

are  vacuum  constants  (p^  = 4 tt 

• 10_;,  eQp0  - 1/c2) 

and  where  terms  0(u2/c2)  are  neglected  since  we  deal  with  nonrela- 
tlvistic  flowfields. 


u « c 


(3) 


In  addition,  we  restrict  consideration  for  the  Maxwell  field  to 
length  and  time  scales,  L and  T,  such  that 


k«  c 

T 


(A) 


the  usual  MHD  limit  that  bars  EM  waves.  Then,  Eq.  la  gives 

E/B  ~ L/T  « c so  that  Eq.  2c  becomes  B = pH  to  0(uL/c2T).  Compar- 

x-w  *>+ 

ing  3D/3t  and  p u to  curl  H,  we  have  in  Eqs.  lc  and  2a  that 

*«sx  6 »v 

curl  H ■=  j to  0(L2/c2T2)  and  j * o(E  + iiu  x K)  to  0(uL/c2T).  In 

xw  /V  XV  /V  XV  /V 

other  words,  we  can  neglect  "flow  magnetization,"  displacement 
current,  and  convection  current  in  Eqs.  1 and  2,  and  they  remain 
valid  up  to  second  order  in  the  smallness  assumptions  of  Eqs.  3 
and  4.  None  of  the  surviving  terms  is  of  first-order  smallness; 
e.g.,  the  "flow  polarization"  in  Eq.  2b  as  well  as  the  induction 
current  in  Eq.  2a  are  each  of  relative  order  u/(L/T). 

A governing  set  of  equations  that  determines  each  CM  field 
variable  individually  as  a functional  of  the  flowfield  u emerges 
from  Eqs.  1 and  2 on  the  assumptions  of  Eqs.  3 and  4.  We  consider 
JB  as  the  primary  field  variable  and  eliminate  IB,  JS,  and  j between 

Eqs.  la,  lc,  2a,  and  2c  to  obtain  the  basic  HM  equation 


3H  . 

~ * — AH  r curl  (u  x H)  , 
3t  po  ~ ~ 


(5a) 


taking  the  ocean  n and  a to  be  constant.  In  addition  to  Eq.  5a,  H 
is  subject  to  the  subsidiary  condition  from  Eqs.  lb  and  2c, 


div  H 


(5b) 
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and  to  the  boundary  conditions  that  follow  from  Maxwell's  Eqs.  lb 
and  lc  (Ref.  12), 

n x H and  n * (F  + M)  continuous  , (6a, 6b) 


where  n is  boundary  normal,  M denotes  magnetization  of  adjoining 
media  (e.g.,  steel  hull),  and  surface  currents  are  assumed  negli- 
gible. Finally,  the  remaining  OM  variables  can  be  written  in  terms 
of  H and  u through  Eqs.  1 through  4: 


B = yH  , (7a) 

rv  a* 

j = curl  H , (7b) 

'»  A/ 

E = — curl  H - yu  x H , (7c) 

/V  0 /N/  A/  - 

D = — curl  H - e0y0  u x H , (7d) 

O A>  A>  Arf 

p = eQyQ(u  • curl  H - H • curl  u)  , (7e) 


where  an  inhomogeneity  j • V(e/o)  was  dropped  in  Eq.  7e. 

The  fundamental  relation  of  the  OM  field  to  the  flowfield 
is  therefore  given  by  Eq.  5a,  which  has  the  form  of  a vector  diffu- 
sion-convection equation.  The  diffusion  term  has  time-scale 
~ ycrL2  and  the  convection  term,  xc  ~ L/u.  Hence, 


convection-rate  d 

diffusion- rate  t 

c 


youL 


*M 


(8) 


where  % is  a magnetic  Reynolds  number.  R^j  is  then  a measure  of 
the  relative  effect  that  the  flowfield  has  on  the  magnetic  field. 

In  most  cases  to  be  considered  in  OM,  Rj^  is  somewhat  less  then 
unity,  though  not  negligible.  Consequently,  the  OM  fields  can  be 
treated  as  small  perturbations  to  the  ambient  field,  though  still 
well  above  the  sensitivity  limits  of  present  magnetometers. 

We  remark  that,  in  the  MHJ)  approximation,  the  familiar  dis- 
placement/conduction current-ratio  tue/o  ~ 10-9/T  is  <<  1/L  by  Eq.  4 
and  so  is  assured  small  for  the  lengths  L considered  here.  That 
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is,  we  do  not  enter  the  electrohydrodynamic  regime  that  obtains 

when  electric  relaxation  becomes  slower  than  magnetic,  viz. 

a/e  ^ 1/t  , or  L ^ 9/uac  ~ 30" 2 meters  here, 
a 


2.2  THE  MAGNETIC  FIELD  ABOVE  THE  OCEAN 

The  utility  of  airborne  observation  warrants  specific  con- 
sideration of  the  field  equations  above  the  sea  surface.  Since 
the  conductivity  of  air  is  ~ lO"14  times  that  of  seawater,  the 
conduction  current  is  completely  negligible  there.  Thus,  within 
the  approximations  of  Eqs.  3 and  4 all  currents  are  zero,  so  that 
curl  H * 0 (cf.  Eq.  7b)  and  we  may  write 


H = grad  V . (9) 


Since  div  H = div  _B/p0  = 0 by  Eq.  lb,  the  potential  y satisfies 
Laplace's  equation 


Ay  = 0 , 


(10a) 


subject  to  the  boundary  conditions  (Eqs.  6a  and  6b),  which  become 
here 


H continuous  at  air/sea  interface 

<v 


(10b) 


and  subject  also  to  the  condition  that  the  field  remain  finite  at 
infinity. 

Therefore,  the  0M  field  above  the  ocean  is  essentially  de- 
termined, via  Laplace's  equation,  by  the  OM  field  at  the  surface 
of  the  ocean.  An  immediate  corollary  is  that  the  altitude  depen- 
dence of  the  OM  field  is  fixed  by  the  horizontal  variation  of  the 
sea  surface  field.  The  generality  of  this  feature  of  the  aero- 
field  is  somewhat  obscured  among  details  in  the  existing  analyses 
cited  in  Sec.  1. 

For  example,  consider  a plane-wave  Fourier  component  of  the 
OM  field  at  the  sea  surface, 


B , - C 

--surface  ~ 


(11a) 
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wmmm 


where  £ is  Fourier  amplitude,  w frequency,  and  horizontal  wave- 
vector.  The  corresponding  Fourier  transform  of  the  aerofield 

(z  > 0)  must  then  have  the  form  K . « V¥  = C e*"wt  ^ ^ z to 

satisfy  Eqs.  10,  wherein  Laplace's  Eq.  10a  requires  A¥  « <2  + q2  =»  0, 
or  q = ±itc.  Hence,  we  have  at  once  that 


H . 
~air 


= H 


-kz 


surface 


(z  > 0) 


(lib) 


taking  Re&  (the  real  part  of  <)  > 0 and  dropping  the  exponentially 
growing  solution.  In  other  words,  a sinusoidal  surface  wave  leads 
to  an  exponential  altitude  decrease  (cf.  Refs.  6 and  8).  More  gen- 
erally, in  a surface  variation  consisting  of  a superposition  of 
Fourier  components,  the  longer  wavelength  effects  tend,  ceteris 
paribus,  to  reach  to  higher  altitudes. 


I 

I 

! 


2.3  THE  OCEAN  FLOW  FIELD 

Returning  to  the  sea,  we  describe  its  flow  field  u by  the 
Navier-Stokes  equation 

du 

p ~~  - -VP  + p(g  - fi  x u)  + nAu  + p(H  • V)H  , (12a) 

wherein  p is  seawater  density,  P is  the  sum  of  kinetic,  magnetic, 
and  viscous  pressures  (P  = p + JpH2  - (5  + n/3)div  u) , g is  earth's 

o> 

gravity  (including  centripetal  acceleration) , ft  its  rotation 

(twice  angular  velocity),  and  n,  C are  shear,  bulk  viscosities 
that  have  been  assumed  constant.  Along  with  Eq.  12a  we  have  the 
continuity  equation 


div  u = 


dlnp 

dt 


(12b) 


and  equations  of  state  and  energy  as  needed. 

The  salient  part  of  Eq.  12a  for  OM  phenomena  is  the  last 
term  together  with  the  magnetic  pressure  term.  These  terms  enter 
via  the  EM  stress  tensor  on  dropping  electric  components 
(E2/c232  ~ L2/c2T*  <<  1 by  Eq.  4)  and  with  p taken  constant  (ambi- 
guities for  moving  magnetic  media  will  not  concern  us  here  — cf. 


- 15  - 


SECTS" 


THE  JOHNS  HOPKINS  UN.vERSlT> 

APPLIED  PHYSICS  LABORATORY 

LAUhEi.  MaPvlmND 


Ref.  13,  p.  496).  Alternatively,  since  V i-H2  = H x curl  H + (H  • V)H, 

»V  /V  »V  »V 

the  two  terms  are  equivalent  to  the  Lorentz  force  density 
p E + j x B « |ij  x H,  to  the  same  approximations  as  Eqs.  5 through 

7.  Notice  that  the  j x H form  makes  clear  that  a strictly  external 

magnetostatic  field  (curl  H « j =0)  contributes  negligible  accel- 

«V  A/ 

eration  to  Eq.  12a.  Hence,  aside  from  any  nonstatic  part  of  the 
ambient  fields  (e.g.,  geomagnetic  time  variations),  magnetism  will 
affect  the  flow  only  through  the  flow- induced  OM  field.  Since  the 
latter  is  expected  to  be  relatively  small  (~R^,  cf.  the  discussion 
following  Eq.  8),  we  therefore  expect  the  effects  of  magnetism  on 
flow  to  be  correspondingly  small. 

Furthermore,  the  relative  size  of  the  magnetic  forces  in 
Eq.  12a  generally  depends  on  the  Alfvdn  speed  A,  which  gauges  the 
magnetic  energy  per  unit  mass.  Thus,  the  ratio  of  magnetic  to 
inertial  terms  in  Eq.  12a  gives 


magnetism 

flow 


where  Az  = 


2 = 


Comparison  of  magnetism  with  other  forces  in  Eq.  12a  can  be  made 
in  terms  of  (A/u)  2,  and  the  usual  Reynolds  number  Re  = uLp/n, 
or  various  combinations.  (For  example,  from  the  last  two  terms  in 
Eq.  12a,  magnetism/viscosity  <*(A/u)^  • Re  or  Rh2/%»  where 

Rjj  5 AlVpop/ri  is  the  Hartmann  number;  we  shall  see  that  the  latter 
measures  the  ratio  of  magnetic  to  viscous  wave  damping.)  For  the 
geomagnetic  field  in  ocean  water,  A ~ 2 • 1CT5  m/s,  which  is  small 
compared  to  most  flow  speeds  of  interest  here.  The  already  small 
OM  magnetism  then  has  its  effect  on  the  flow  reduced  by  the  small 
ratio  (A/u)2  of  Eq.  13. 

Consequently,  it  seems  well  justified  to  ignore  the  magnetic 
fc rces  in  Eq.  12a.  This  neglect  uncouples  Eqs.  12  from  Eqs.  5,  and 
ths  flow  description  becomes  an  ordinary  hydrodynamic  problem  that 
we  shall  assume  to  have  been  solved  elsewhere.  We  are  then  left 
to  consider  Eqs.  5 alone,  as  outlined  in  Sec.  3.1  and  elaborated 
in  the  bulk  of  this  report  (Secs.  4 through  6).  Because  of  the 
potential  importance  of  reciprocal  magnet ism- flow  interactions  as 
propagators  of  OM  waves,  we  shall  also  look  more  closely  at  the 
coupled  Eqs.  12  and  5 in  Sec.  3.2  under  various  simplifying  approxi- 
mations. 


Ref.  13.  R.  M.  Fano,  L.  J.  Chu,  and  R.  B.  Adler,  Electro- 
magnetic Fields,  Energy  and  Forces,  Wiley,  New  York,  NY,  1960. 
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Table  1 presents  some  definitions  and  numerical  values  of 
various  parameters  of  importance  for  ocean  magnetic  phenomena. 


Table  1 

Ocean  Magnetics  Parameters 


Symbol 

Parameter 

Nominal  Value* 

a 

Ocean  electrical  conductivity 

4 

V 

Ocean  magnetic  permeability 

4ir  x 10“ ' 7 

t 

Ocean  dielectric  constant 

9 x 10-9/4tt 

6 

Ocean  compressibility 

4.5  x 10"10 

n 

Ocean  shear  viscosity 

1.5  x 10“3 

C 

S 

Ocean  sound  speed 

1500 

p 

Ocean  density 

1025 

H 

e 

Earth  magnetic  field 

5 x 104  (Y) 

8 

Earth  gravity 

9.8 

*mks  units  are  used  throughout,  except  that  magnetic  fields  are 
measured  in  y (=  1C-5  G). 
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3.  OVERVIEW  OF  PHENOMENA 


We  decompose  the  magnetic  field  and  the  flow  field  into 
their  uncoupled  portions  plus  hydromagnetically  coupled  remainders. 


H = Kn  + h , u = uft  + v . 

/V  AW  AW  Aw*  AW 


(14a, 14b) 


That  is,  represents  the  magnetic  field  neglecting  flow  and 
the  flowfield  neglecting  magnetism,  both  of  which  will  ordinarily 
be  assumed  as  given  here. 

The  OM  field  proper,  h,  v,  should  be  small  according  to  the 
discussions  following  Eqs.  8 and  13;  specifically,  we  expect  the 
hierarchy 


-*■«£«  1 . 
«o  Ho 


(Or,  if  »£o  m 0,  consider  the  first  ratio  to  be  v/(L/T)).  There- 
fore, in~the  first  instance  one  is  led  to  consider  the  limiting 


v = 0 (OM  induction)  , 


(16a) 


in  which  & is  simply  induced  by  the  given,  purely  mechanical,  flow 
Uq.  Subsequently  we  may  lift  this  restriction  in  order  to  investi- 
gate the  important  (however  unlikely)  possibility  of  significant 
OM  wave  propagation.  That  is,  the  general  case 


v / 0 (OM  radiation) 


(16b) 


admits  h and  v interactions  that  may  conceivably  radiate  through 
the  ocean. 
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3.1  VARIETIES  OF  OM  INDUCTION 

i In  the  case  given  by  Eq.  16a,  the  OM  field  h is  essentially 

i determined  via  the  basic  HM  Eq.  5a  as  a function  of  the  given 

fields  H0,  Uq.  rhis  gives  a straightforward  induction  problem, 
which  becomes  multifarious,  however,  through  the  variety  of  non- 
! hydromagnetic  fields  Ho  and  uo  that  will  enter. 

j *** 

The  given  magnetic  field  is  taken  as  tne  sum  of  earth  (e) 

| and  magnetized-body  (b)  fields, 

, 50-2e+j^  * (17) 


Hg  being  the  total  geomagnetic  field  and  Hj>  the  proper  field  of 
the  body  when  immersed  at  rest  in  quiescent  seawater. 

In  general,  Hg  consists  of  the  magnetostatic  mean  earth 
field,  which  can  be  taken  as  locally  uniform  over  scales  of  in- 
terest here,  plus  the  slight  temporal  and  local  variations  often 
called  geomagnetic  and  geologic  noise,  respectively.  (Notice  that 
the  noise  includes  that  generated  by  currents  induced  in  the  sea 
by  the  geomagnetic  noise  itself,  so  that  He  is  not  a strictly  ex- 
ternal field,  cf.  Eq.  20a  below.)  The  noise  is  usually  well  above 
magnetometer  sensitivities  and  therefore  must  be  considered  in  the 
detection  of  signals  such  as  the  OM  h field  down  to  the  limits  of 
sensitivity.  However,  the  hydromagnetic  effects  of  the  noise  are 
of  course  much  smaller  than  those  of  the  mean  geomagnetic  field. 
Therefore,  we  can  set  aside  the  geomagnetic  and  geologic  noise  in 
the  calculation  of  the  OM  £ signals,  and  consider  Hg  as  constant 
in  space  and  time. 

Similarly,  while  in  general  has  a spatial  variation  that 
we  shall,  of  course,  not  ignore,  its  intrinsic  time  variations 
will  be  assumed  minor  here.  (Notice  that,  like  Hg  above,  Hj,  in- 
cludes a part  associated  with  currents  induced  in  the  seawater  by 
any  intrinsic  time  variation  of  ^ itself,  so  H^  too  is  not  a 
strictly  external  field,  cf.  Eq.  20b  below.)  On  that  assumption, 
we  can  neglect  the  hydromagnetic  effects  of  the  intrinsic  time 
variations  in  calculating  ji  and  consider  Hj,  as  constant  in  the 
body  rest  frame.  There  remains,  of  course,  its  convective  time- 
variation  due  to  motion  of  the  body  through  the  sea  (cf.  the  last 
term  in  Eq.  22a). 

The  given  mechanical  field  is  taken  to  be  composed  as 

u0  = "natural"  and/or  "body-related"  motions  , (18) 
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wherein  the  two  classes  of  motion  include 
1.  The  natural  motions. 

Surface  gravity  waves  (SW) , due  to  either  local 


winds  (sea)  or  distant  storms  (swell) , (18a) 
Internal  gravity  waves  (IW) , (18b) 
Ambient  acoustic  waves  (AW),  and  (18c) 
Sea  currents  and  tides;  (18d) 


2.  The  body-related  motions  due  to  passage  of  a vehicle 

through  the  seawater: 

Mean  flow,  i.e.,  the  ~elative  translation  of 

the  seawater  and  the  body-field  H^,  (18e) 

Potential  flow,  caused  by  hull  displacement,  (18f) 

Wake  flow,  from  turbulent  growth  near  the  body 

to  collapse  in  far  region,  and  (18g) 

Propagated  disturbances,  due  to  body-generated 

SW,  IW,  and  AW.  (18h) 

In  general,  these  various  motions  are  not  independent,  due  to  hy- 
drodynamic interactions,  so  that  Eq.  18  is  a nonadditive  combi- 
nation. However,  for  simplicity  we  will  consider  the  idealiaed 
situation  of  negligible  interaction.  Then  Eq.  18  becomes  additive, 
as  does  Eq.  17,  so  that  the  hydromagnetic  driving  (last  term  of 
Eq.  5a)  is  a sum  of  binary  combinations  of  each  part  of  H0  with 
each  of  the  uQ. 

The  numerous  possible  pairs  of  given  motion  (Eqs.  18a  through 
18h)  and  magnetism  (Eq.  17)  yield  a respectably  long  list  of  phenom- 
ena to  be  considered.  After  a general  theoretical  development  in 
Sec.  4,  each  phenomenon  will  be  considered  individually  in  Secs.  5 
and  6.  But  first  we  will  discuss  briefly  the  propagation  of  purely 
inductive  effects  (pseudowaves)  versus  true  OM  radiation. 


3.2  PSEUDOWAVES  AND  OM  RADIATION 

Even  in  the  case  v = 0 (Eq.  16a)  of  pure  Oh  induction 
notice  that  there  can  be~propagation  of  the  induced  magnetism  h 
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when  the  inducing  flow  is  a mechanical  radiation.  That  is,  if  the 
prescribed  uQ  is  a propagating  flow  (SW,  IW,  or  AW),  then  the 
purely  mechanical  waves  carry  along  with  them  their  induced  h 
fields.  Such  OM  "pseudoradiations"  are  treated  in  detail  as 
propagating  varieties  of  pure  induction  phenomena  (Eqs.  18a,  b,  c, 
and  h,  above)  in  Secs.  4 through  6. 

In  the  case  of  v j*  0 (Eq.  16b),  however,  one  has  the  possi- 
bility of  the  complete  OM  fields,  £ and  v,  mutually  interacting  so 
as  to  produce  a magnetomechanical  wave  propagation,  or  true  OM  ra- 
diation. In  that  case,  we  have  to  consider  the  coupled  set  of  Max 
well-Navier-Stokes  equations  that  we  reduced  in  Sec.  2 essentially 
to  Eqs.  5 and  12  plus  subsidiary  relations  as  required.  These 
equations  become  by  virtue  of  Eqs.  14  and  15  a linearized  set  cou- 
pling h and  v.  With  h considered  determined  in  terms  of  v by  Eqs. 

5 and  with  g fl  specified,  Eqs.  12  represent  four  relations  among 
~»  ~ 

five  unknowns  (v  and  the  OM  perturbations  in  density  and  pressure, 
6p  and  6p) . To  complete  the  system,  simple  equations  of  state 
and/or  energy  will  suffice  for  the  idealized  situations  (cases  1, 
2,  and  3,  below)  considered  here.  More  general  cases,  of  course, 
would  require  more  detailed  considerations  of  state  and  energy 
conservation. 


1.  Incompressible  and  Homogeneous  Ocean:  Under  the 
simplest  assumption  the  state  equation  is 


p = Constant  in  Space  and  Time 

In  this  case,  Eqs.  5 and  12  yield  the  familiar  Alfven  radiation 
(e.g.,  Ref.  14).  For  example,  neglecting  viscosities  and  earth 
rotation  ii.  Eq.  12a  (Eq.  12b  is  simple  div  _u  = 0)  and  supposing 
the  unperturbed  ocean  is  at  rest,  we  get  the  plane-wave  dispersion 
relationship 


= k2A2cos20jcy 

between  wave-vector  (k)  and  frequency  (w) , with  A the  Alfven  speed 
(Eq.  13)  and  0^  the  angle  between  k and  Hq.  But  in  the  ocean,  7 
and  A are  so  small  that  the  dissipative  term  predominates,  and 


f ik2l 

> b> 

l wo 


Ref.  14.  J.  M.  Robertson,  Hydrodynamics  in  Theory  and  Appli- 
cation, Prentice-Hall,  NJ,  1965. 
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k2  w -ipoto  at  any  sensible  frequency  (viz.,  to  >>  poA2  ~ 10-11  Hz). 
In  other  words,  an  Alfv^n  wave  is  so  highly  damped  in  the  ocean 
that  it  degene;,  tes  to  simply  a magnetic-diffusion  wave  (solution 
of  Eq.  5a  with  u ■ 0). 


2.  Incompressible  but  Inhomogeneous  Ocean:  If  the  ocean 
is  treated  still  as  incompressible,  but  inhomogeneous,  we  can  con- 
sider the  density  a Lagrangian  constant,  i.e.,  the  equation  of 
state  is 


ie. 

dt 


or 


= -(u 
9t  ~ 


V)p 


In  this  case,  Eq.  12b  is  again  div  u * 0,  but  Eq.  12a  with  vis- 
cosities and  earth  rotation  neglected  involves  the  added  term  g6p, 

leading  to  internal  gravity-magnetism  waves.  For  example,  if  the 
unperturbed  state  is  one  of  stable  vertical  (z)  stratification, 
then  the  above  equation  linearizes  to  (9/3t  + u0  • V)6p  = -v2  dp0/dz. 

For  a constant  buoyancy  frequency  N (i.e.,  p q “ exp(-N2L/g))  and  in 
the  Boussinesq  approximation,  one  obtains  the  plane-wave  dispersion 
relationship 


a2  - N2sin20 


id  - 


ik2l 

pa 


(DA2k2cos20 


kH 


where  0 is  the  angle  k makes  with  the  vertical.  Here  if  one  ne- 
glects the  oceanic  A and  a as  vanishingly  small  (again  for  fre- 
quencies >>  10-11  Hz),  one  obtains  in  addition  to  the  degenerate 
Alfv£n  (magnetic-diffusion)  wave  above,  the  usual  internal-gravity 
modes  (id  » ± N sin0)  for  constant  N.  A more  complete  discussion, 
which  also  includes  shear  flow  effects,  will  be  found  in  Ref.  15. 


3.  Compressible  but  Adiabatic  Ocean:  When  one  admits 
compressible,  but  adiabatic,  ocean  motions,  the  linearized  equa- 
tion of  state  is 


p = p0(l  + 0(p  - p0)]  or  5p  = 


cs2 


Ref.  15.  N.  Rudraiah  and  M.  Venkatachalappa,  "Effect  of 
ohmi;  dissipation  on  internal  AlfvSn- gravity  waves  in  a conducting 
shea.:  flow,"  J.  Fluid  Mech. , Vol.  62,  1974,  p.  705. 
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where  6 = l/pQCg2  is  the  compressibility  and  cg  the  sound  speed  in 

the  sea.  In  this  case  Eqs.  5 and  12  yield,  in  addition  to  the 
degenerate  Alfven  wave  of  case  1 and  the  IW  for  stratified  in- 
homogeneity of  case  2,  the  magnetosonic  radiations  of  compressible 
MHD.  Neglecting  damping  and  assuming  pQ  to  be  uniform,  we  have 
for  the  dispersion  relation  of  the  magnetosonic  modes 


■i  k2{<% 


-)  + 


(c  2 + A2)2  - 4c  2A2cos20, 

S S Kn 


But  in  the  ocean  A2/cs2  ~ 10“12,  so  that  with  very  little  error 
the  magnetosonic  wave  is  simply  the  ordinary  acoustic  compression 
wave  (id2  = k2cs2).  More  detailed  considerations,  which  also  show 
the  very  slight  MHD  effects  on  damping  and  acoustic  shear  waves, 
will  be  found  in  Ref.  16. 

From  these  examples,  it  will  be  clear  that  true  OM  radia- 
tion is  not  expected  to  be  of  significance  in  the  ocean.  However, 
we  emphasize  that  the  purely  mechanical  radiations  (SW,  IW,  and 
AW)  to  which  the  magnetogravity  and  magnetosonic  modes  reduce  in 
the  ocean  are  still  capable  of  propagating  induced  magnetic  fields, 
the  OM  pseudoradiation  referred  tc  above.  Thus,  for  example, 
ocean  compressibility  will  allow  what  we  will  call,  in  contradis- 
tinction to  the  magnetosonic  waves  above,  the  sonomagnetic  OM 
pseudowaves  discussed  in  Subsecs.  5.3  and  6.6. 


Ref.  16.  G.  I.  Cohn,  "MHD  wave  phenomena  in  seawater,"  Elec- 
tromagnetics of  the  Sea,  AGARD  Conf.  Proc.  77,  1970. 
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4.  THEORY  OF  OM  INDUCTION 
AND  USEFUL  FORMULAE 


In  the  pure  induction  case  (Eq.  16a)  in  which  the  flow  is 
prescribed  (u  = u0),  Eqs.  5 alone  determine  the  magnetic  field, 


H = H.  + h = H + HL  + h (19) 

'V  *v»V  ««w  /V 

(cf.  Eqs.  14a  and  17).  By  our  definition  of  Hq  as  independent  of 
flow,  we  have  and  satisfying  Eqs.  5 in  the  forms 

3H 


AH 

!% 

o 

TX 

II 

(20a) 

fa  1 

= vah  = V0{K+Zb  *IJ»b  * 

(20b) 

div  H 
~e 

= div  = 0 . 

(20c) 

In  Eq.  20b,  ^ is  time  derivative  in  the  body  rest  frame,  i.e.,^ 
is  to  be  independent  also  of  the  mean  relative  flow  of  seawater 
corresponding  to 


arbitrary  body  velocity  = V.  (t) 


(21) 


One  sees  that  in  general  neither  Hg  nor  Hb  are  external  fields, 
i.e.,  A Hp  b 5^  0,  since  their  time  variations  induce  currents  in  the 
ocean  ^ r.  curl  ^ ^ ^ 0.  As  indicated  in  Sec.  3,  we  simply  as- 
sume the  solutions  ^ of  Eqs.  20  as  given  for  present  purposes. 

On  inserting  the  decomposition  (Eq.  19)  for  H into  Eqs.  5 
and  using  Eqs,  20,  we  obtain  the  determining  equations  for  the  OM 
field, 


f 

at 


po 


h = curl(ur  x H) 


<S> 


(22a) 


div  h = 0 


(22b) 


which  are  exact  within  the  HM  approximation  (Eq.  4)  and  the  induc- 
tion assumption  (Eq.  16a). 


i 


! 

■n-4r 
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An  approximate  form  of  Eq.  22a  will  suffice  for  discussion 
here.  The  convective  term  curl  Qjq  x $)  is  expected  to  be  small 
(cf.  Eq.  8),  so  that  we  shall  set  jj  « HQ  in  it.  Then,  symmetrizing 
Eq.  22a  by  vector  identities  and  Eq.  20c,  we  have 

(a  - pa  -^)h  « yo  curKHoX  u0  - HfoX  . (23) 

Since  this  is  linear  in  we  can  handle  the  body  magnetization- 

translation  interaction  x V^,)  separately  from  the  other  inter- 
actions (Hq  x u0)  below.  In  the  latter,  we  usually  consider  HgwHe 
(i. e. , <<  /SeT  Tn  line  with  the  remarks  following  Eq.  17,  we 

treat  ^ as  uniform  and  static  and  Hj,  as  intrinsically  static, 
while  of  course  Vj>  is  spatially  invariant.  Therefore,  when  V^,  is 
steady,  both  HQ  « Hg  and  Vjj  can  be  considered  absolute  constants. 
Notice  that  the  flow  ^ is  unrestricted  (e.g.,  div  ^ / 0)  and 
that  h <<  Ht.  is  not  assumed. 

/v  *wD 

In  general,  to  determine  the  complete  march  of  magnitude 
and  configuration  for  the  vector  field  h requires  a detailed  analy- 
sis for  each  of  the  many  varieties  of  given  flow  uq  (Sec.  3).  Var- 
ious approaches  are  appropriate  for  various  flows'?  as  indicated 
briefly  below  (and  in  Appendix  A).  However,  before  carrying  out 
detailed  analyses,  one  would  like  to  have  some  estimates,  however 
crude,  for  the  0M  fields.  Estimates  for  the  sizes,  and  somewhat 
for  the  extents,  are  readily  derived  and  will  be  the  main  subject 
of  the  rest  of  this  report. 


4.1  LOCAL  ESTIMATOR  FOR  THE  SIZE  OF  OM  FIELDS 

A rule  of  thumb  for  the  magnitude  of  the  OM  field  induced 
locally  follows  from  a dimensional  analysis.  Thus,  setting 
9/3t  ~ 1/T,  A2  ~ 1/L2,  and  curl  ~ 1/L  (cf,  Eq.  4)  in  Eq.  23  with 
the  body  term  put  aside,  we  get 


UOu0LH0 

h ~ 1 + 0(id/T)  for  T » td  ’ *M  ‘ H0  (24 

(cf.  Eq.  8).  Here  L is  the  flow  scale,  which  is  assumed  to  set 
the  h field  scale.  (For  the  body  term  in  Eq.  23,  let  uQH0 

and  L be  the  Hb  scale.)  The  last  form  in  Eq.  24  is  the  quasi- 
static limit  (rate  of  change  1/T  « diffusion  rate  1/Td),  which 
will  be  applicable  for  most  OM  induction  effects.  The  opposite 
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limit,  T <<  t,  or  L2  >:>  T/pcr  for  which  h/Ho  ~ R^T/t^  = uq/(L/T), 
becomes  important  only  for  very  long  wave  motions  (cf.  Sec.  4.3, 
where  one  notes  that  "quasi-static"  does  not  mean  "low  frequency" 
for  oceanic  waves). 

Numerically,  the  estimator  (Eq.  24),  with  Hq  « ^ and  the 
values  in  Table  1,  reads 


1 -5 
h ^ UgL  (in  y = 10  G) 


for  uq  in  m/s  and  L in  m.  We  emphasize  that  this  is  only  an  order- 
of-magnitude  upper  limit,  and  in  special  cases  the  field  could  be 
much  smaller.  (For  example,  for  a uniform  field  jHo  and  an  incom- 
pressible flow  juq  that  varies  only  transversely  to  Hq,  the  right- 
hand  side  of  Eq.  23  vanishes  and  ih  a 0.)  Nevertheless,  Eq.  25 
affords  a useful  rule  of  thumb  when  used  with  proper  discretion. 


4.2  A THEOREM  ON  THE  EXTENT  OF  OM  FIELDS 

For  a localized,  incompressible  flow  in  the  quasi-static 
limit,  an  explicit  formula  for  the  OM  field  h at  distant  points  K 
can  be  derived  in  terms  of  the  spin  angular  momentum  £ of  the  flow, 
as  follows. 

Equations  23  and  22b,  with  9/3t  neglected  (T  >>  tj),  H0s*,ge 
uniform,  div  uQ  = 0,  and  the  body  term  ret  aside,  possess  the  inte- 
gral solution 


h(R)  = j-2.  (H 
~ ~ 4 it  ~e 


V)  f dr  uQ 

/■w  «r  «vv 


plus  vector  harmonics  as  needed  to  satisfy  boundary  conditions 
(cf.  Appendix  A).  Consider  distant  field  points  Swell  outside 
the  flow  region.  Then  the  integral  in  Eq.  26  can'Tie  approxi- 
mated by  a Maclaurin  expansion  of  1/1 R - r!  as 

fdz  u0i7ls  " = (1/R)/  dl  uoi  - S dl  rj  uoi  + oa/R3) 

in  tensor  notation.  In  these  integrals,  by  hypothesis: 
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1.  div  uQ  = 0,  so  that  uQi  = div(r^  uQ)  and  tj  u0^  = 
IdivC^r^  u0)  + |(uoirj  - u^r.^,  by  vector  identities; 
and 

2.  u is  localized,  so  that  the  divergence  terms  give  van- 
ishing surface  integrals  (specifically,  we  require 

to  decrease  faster  than  l/r1*  asymptotically). 

Thus,  we  get 

f dr  u0/|R  - r | » \ (V  1/R)  x fdr  (r  x u0)  , (27) 

in  which  we  recognize  the  last  integral  to  be  S/p  for  a flow  re- 
gion of  total  spin  £ and  homogeneous  density  p.  Combining  Eqs.  27 
and  26  gives  then 

h(R)  = - ~ (K  * f)(S  x V)  j+O  k (28) 

for  sufficiently  localized,  incompressible,  quasi-static  flow  in 
approximately  homogeneous  infinite  medium.  The  presence  of  the 
air-sea  boundary  will  modify  this,  but  the  complete  solution  will 
retain  comparable  magnitude  in  general.  (An  exception  is  the 
highly  symmetric  case  where  ,§  and  jgg  are  both  vertical,  for  which 
the  complete  solution  includes  an  image  field  that  exactly  cancels 
Eq.  28  everywhere  on  the  sea  surface  anc',  thus,  bv  Eq.  10b  gives 
zero  aerofield;  cf.  Sec.  6.4.) 


1 


Therefore,  on  working  out  the  derivatives  in  Eq.  27,  we 
have  the  explicit  theorem 


h(R)  «f-Sx[H  - 3R(R 


£*»  h 


under  the  conditions  stated  above,  wherein  jl  = R/R.  Notice  that 
this  falls  off  with  distance  like  a dipole  field,  « 1/R3,  but  of 
course  with  a more  complicated  directional  dependence.  Its  magni- 
tude Mien  may  be  represented  in  terms  of  effective  dipole  magni- 
tude, 
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where 


m 


eff 


VioSHe 

8irp  ‘*e 


Vol 

8ti 


(30b) 


in  which  Vol  denotes  the  flow  region  volume.  Alternatively,  one 
can  write  the  distant  field  estimate  (Eqs.  30)  as 


^(distant)  ^(local) 


' Vol ' 
8ttr3 


(31) 


in  terms  of  the  local  estimator  of  Eqs.  24  and  25. 

The  above  theorem  and  its  derivation  are  mathematically 
analogous  to  a well-known  magnetostatic  theorem  (Ref.  17,  p.  145) 
expressing  the  vector  potential  due  to  localized,  divergenceless 
currents  in  terms  of  their  magnetic  moment:  A 11  sJ  - 

(V.  1/R)  x m (cf.  Eq.  27).  However,  there  is  the  important  dif- 
ference that  the  induced  OM  currents  here  are  not  localized.  That 
is,  whereas  in  the  magnetostatic  case  the  exterior  current  is  zero 
(g_  « curl  A « grad  (m  • V)l/R,  so  that  j = curl  H = 0) , in  the  in- 
duction case  here  one  has  from  Eqs.  7a  and  24  the  distant  (’’return”) 
current 


j(R)  « 


HO 

37Tp 


<5. 


V)(S 


(32) 


which  is  nonzero. 

The  theorem  29  and  resulting  estimates  30  through  32  will 
clearly  be  most  useful  for  localized  flows  such  as  an  isolated 
wave  crest  or  the  wake  of  a body  traversing  the  ocean. 


Ref.  17.  J.  D.  Jackson,  Classical  Electrodynamics,  2nd  edi- 
tion, Wiley,  New  York,  NY,  1975. 
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4.3  FLOW-WAVE  INDUCTION  OF  OM  PSEUDOWAVES 

At  the  other  extreme  from  localized  flows,  consider  flow- 
fields  of  large  spatial  extent  that  can  be  Fourier  space-time 
transformed  into  ^armonic  wave  components.  In  the  ocean,  familiar 
examples  are  various  gravity  waves  and  acoustic  fields.  Then,  one 
has  to  deal  with  induction  by  the  transformed  flow  in  Eq.  23  that 
becomes  linear  with  constant  coefficients  on  setting  aside  the 
body  term  x V^,  taking  Hg  w He  uniform  and  static,  and  pa  con- 
stant. Thus,  Eq.  23  trar  forms  to  yield  the  OM  field  also  in  har- 
monic waves: 

~ . . imt-iic*r  „ , *>•,  , imt-iKT 

u„  ->  un(z)e  ~ ~ =>  h -*•  h(z)e  (33) 

*S#V  /S^  U ~ (Nrf 


for  a flow  wave  of  (angular)  frequency  o>  and  horizontal  wave  vec- 
tor £.  This  assumes  an  ocean  of  infinite,  homogeneous  horizontal 
extent,  but  of  course  vertical  (z)  discontinuity  or  inhomogeneity 
is  allowed  for  in  the  amplitudes,  ug(z)  and  jj(z).  Notice  that  the 
OM  frequency  and  wave  vector,  and  correspondingly  its  phase  and 
group  velocities,  are  identical  with  those  of  the  flow  wave.  That 
is,  the  induction  wave  propagates  exactly  as  does  the  given  flow 
wave.  Physically,  the  magnetism  simply  rides  along  on  the  mechani- 
cal wave,  whence  we  call  it  OM  pseudoradiation. 


The  amplitude  of  the  OM  pseudowave  (Eq.  33)  may  be  estimated 
from  Eq.  23,  just  as  was  done  Jn  Eq.  24  and  the  following  discussion, 
with  T -*■  l/u>.  Thus  we  have  the  limiting  behaviors 


V He 


— • H 
(oL  e 


wL2  <<  — 


ioL2  » — 


(34a) 


(34b) 


where  R^  = pougL.  Except  for  waves  long  enough  for  the  vertical 
scale  length  to  enter,  we  have  L ~ 1/k  and  wL  ~ the  wave  phase 
speed,  Cph  £ w/k  (cf.  L/T  -*■  Cph  below  Eq.  24).  Equation  34a  corre- 
sponds to  the  quasi-static  limit  in  Eq.  24,  and  Eq.  34b  to  the 
nonstationary  limit  indicated  following  Eq.  24.  Notice  that  in 
the  latter  case  the  OM  field  is  independent  of  the  ocean  param- 
eters p and  o,  except  as  they  delimit  the  to-L  regime.  We  empha- 
size that  the  quasi-static  limit  is  not  generally  the  same  as  a 
low-frequency  limit.  Rather,  for  oceanic  varieties  of  wave  motion. 
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Eq.  34a  is  the  short  wavelength  (X  = 2tt/ic),  high-frequency  limit, 
while  Eq.  34b  obtains  for  long,  low-frequency  waves. 

To  be  specific,  the  two  regimes  in  Eqs.  34  are  delimited  by 
the  criterion  function  to  = 1/poL2,  which  is  tc2/po,  whereas  all 
ocean  waves  have  the  dispersion  relation  u>  = m(j&)  increasing  less 
rapidly  than  <2  for  all  jc.  Consequently,  the  quasi-static  case  (Eq. 
34a)  is  certain  to  obtain  for  u>,  k sufficiently  above  (or  X below) 
a point  mo,  <o  (or  Aq).  On  the  other  hand,  the  opposite  limit 
(Eq.  34b)  may  or  may  not  obtain  at  to  < tog,  X > X0,  depending  on 
the  vertical  scale  length  (e.g. , see  Subsec.  5.1.5.  In  particu- 
lar, the  dispersion  relations  for  SW  on  ocean  of  depth  D,  for  IW 
on  a relatively  sharp  density  change  6p/p  at  depth  d <<  D,  and  for 
AW  of  constant  speed  cs  with  a hard  ocean  bottom,  are  respectively 


SW: 

l>2(k) 

= g«  tanh  kD  , 

(35a) 

IW: 

ui2(k) 

g«(6p/p) 

(35b) 

~ 

[ coth  ted  + coth  k(D 

- d)]  ’ 

AW: 

V 

CM 

3 

» c 2 

(£  + h\Zn2 

k2  + i J 

5 

- 

(modes  l - 0,1,2,...)  . (35c) 


(For  the  shorter,  capillary  SW  and  the  small-scale  IW,  as  well  as 
Coriolis  effects,  see  Ref.  18.  The  hypothetical  case  of  "bound" 
SW  is  discussed  in  Subsec.  5.1.5.)  In  the  three  cases  (Eqs.  35), 
one  finds  for  the  quasi-static  limit  above  (to  > Uq  and  X < X0  = 


2ir//pa(D3  ) the  values: 

- 

SW:  = pogD 

for 

D < 1 km  , 

i 

toQ  = (pog-J3 

for 

D £ 2 km  ; 

IW:  u>0  = pagd  (6p/p)  ; 

AW:  a ).  = pac  2 
0 s 


The  frequency-wavelength  (v  2 u)/2u,  X = 2 it/ k)  dispersion 
curves  (Eqs.  35)  are  illustrated  numerically  in  Fig.  1 for  tvpical 


Ref.  18.  0.  M.  Phillips,  The  Dynamics  of  the  Upper  Ocean, 

Cambridge  University  Press,  1966. 
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Fig.  1 Ocean  Wave  Dispersion  Curves  of  Acoustic  (AW),  Surface  (SW),  and 
Internal  (IW)  Wave  Motions  for  Various  Ocean  Depths  (D),  Pycnocline 
Depths  (d),  and  Density  Changes  (5 p/p),  Sonic  Mode  (£),  with  IW 
Truncated  at  Brunt  Vaisala  Frequency  for  a 25-m-Thick  Pycnocline. 
Dashed  line  separates  quasi-static  (left)  and  nonstatic  (right)  regimes  of 
OM  pseudoradiation  (cf.  Eq.  34a). 
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ocean  parameters.  The  dashed  diagonal  .Hne  represents  the  limiting 
locus  (vg,Ao)  for  the  quasi-static  approximation  (Eq.  34a)  to  be- 
come valid  for  infinite  vertical  scale  length,  while  departures 
from  this  locus  due  to  finite  vertical  scale  are  indicated  at 
k < 1/D  (X  > 2irD) . One  sees  that  the  quasi-static  approximation 
Eq.  34a  holds  for:  (a)  all  internal  waves  but  those  of  tidal  fre- 
quencies; (b)  all  surface  waves  to  well  below  the  swell  band, 
breaking  down  only  for  the  rare  tsunami  or  longer  waves;  and 
(c)  all  acoustic  waves  of  audio  frequency  or  higher.  The  non- 
static limit  (Eq.  34b)  can  obtain,  on  the  other  hand,  for  deep- 
water tsunamis,  tides,  and  infrasound. 
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5.  ESTIMATES  OF  OM  INDUCTION 
BY  NATURAL  SEA  MOTIONS 


From  the  foregoing  theory,  we  will  estimate  in  this  section 
the  OM  fields  induced  by  the  various  flows  due  to  natural  sea  move- 
ments, as  listed  in  Eqs.  18a  through  18d.  Aside  from  their  in- 
trinsic interest  (e.g.,  for  ocean  wave  sf'dies) , these  naturally 
induced  fields  are  important  as  noise  background  for  observations 
on  the  body- related  OM  fields  treated  in  Section  6. 

As  indicated  in  the  Introduction,  the  estimates  in  this  sec- 
tion and  in  Sec.  6 are  intended  to  be  simple  and  clear,  rather 
than  precise.  When  available,  more  detailed  analyses  will  be 
cited  for  more  exact  results. 


5.1  SURFACE  WAVES 


The  OM  field  induced  by  the  interaction  of  Hg  with  the  most 
common  ocean  waves,  sea  and  swell,  have  been  discussed  by  a number 
of  authors  (cf.  Refs.  4,  6,  and  8,  and  others  cited  therein). 
Longer  waves  such  as  tsunamis  were  considered  by  Larsen  (Ref.  9). 
We  can  collate  and  summarize  the  essential  results  of  these  de- 
tailed calculations  by  means  of  the  foregoing  theory.  First  we 
consider  the  cases  of  sea  and  swell  on  deep  ocean,  as  well  as 
swell  in  shelf  waters,  which  come  under  the  quasi-static  approxi- 
mation of  Eq.  34a.  We  then  go  on  to  the  longer  waves  for  which 
Eq.  34b  obtains. 

Linear  plane  gravity  waves  of  frequency  u>  = 2irv,  wave  vec- 
tor £,  and  amplitude  av  have  fJow  speed  and  scale  roughly 


v 


tanh  kD 


and 


L 


tanh  kD 
k 


(36) 


near  the  sea  surface  (z  = 0) , where  of  course  uQ  is  largest.  Then 
one  has 


h ^ H yoa  — 
e v k 


1 ~ 0) 
T a — 
4 V K 


y 


(37) 


from  the  quasi-static  Eq.  34a,  which  applies  to  both  swell  and  sea 
(Fig.  1). 
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On  the  deep  ocean  (D  > 1/ic,  w2  « gx),  "uj  decreases  « e 
below  the  surface  (z  < 0),  and  completing  the  OM  solution  to  fit 
the  boundary  condition  of  Eq.  10b  introduces  a factor  of  5 into 
the  surface  field.  Then  more  accurately 


h 


surface 


(38) 


while  h decreases  essentially  exponentially  with  depth  (~  2xze 
for  -z  >>  1/x).  At  altitudes  z i 0,  Eq.  11b  with  x =»  u>2/g  «*  4v2 
gives  at  once,  from  Eq.  38, 


h . « 

air 


-f-1 

10  M 


e ^ZV  y (deep  water) 


(39a) 


t 

Detailed  analysis  modifies  these  expressions  by  a directional  fac- 
tor involving  the  magnetic  dip  angle  (I)  and  the  wave  propagation 
angle  (0^)  relative  to  magnetic  north,  namely. 


V sin2I  + cos2I  cos20 

X 


y 


which  is  of  order  unity,  except  for  east-west  waves  in  the  tropics 
when  it  approaches  zero  (cf.  Ref.  6).  We  set  aside  that  case  here. 

In  shallow  waters  (xD  < 1,  10  » xi'gD),  bottom  effects  become 
important.  The  flow  amplitude  'u0  has  little  depth  variation,  and 
the  OM  field  penetrates  to  the  seabed  where  it  must  be  made  con- 
tinuous into  the  conducting  sediments.  Detailed  solutions  are 
given  by  Larsen  (Ref.  9)  and  numerical  results  by  Klein  et  al. 

(Ref.  4).  But  observe  simply  that  the  above  estimates  remain 
valid  as  to  order  of  magnitude  provided  that  instead  of  x w w2/g 
we  use  x « u/ /gD,  according  to  Eq.  35a.  Thus,  Eq.  39a  becomes 


*av  e Y (shallow  water) 


(39b) 


1 9 

Note  from  Eq.  35a  that  xD  < 1 implies  D < -r T , which,  for  periods 
T up  to  the  longest  (20  s)  swells,  gives  D 4 80  m,  i.e.,  the  shal- 
lowest shelf  waters. 


The  above  estimates  make  clear  the  very  few  factors  tnat 
fix  the  gross  field  magnitude.  The  dominant  factor  is  the  exponen- 
tial altitude  (and  depth)  behavior,  such  that  the  lowest  frequency 


S 1 


| 
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waves  can  be  larger,  since  ay  « 1/v2  for  given  steepness  (2av/A). 
Thus,  the  lowest  frequency  waves  also  tend  to  induce  OM  fields  of 
the  most  sizable  amplitude.  More  precise  statements  require  speci- 
fication of  the  wave  power  spectrum,  which  depends  on  local  winds 
and  distant  storms.  Representative  spectra  fcr  the  corresponding 
sea  waves  and  ocean  swells  are  sketched  in  Fig.  2.  We  stress, 
however,  that  the  actual  spectra  obtaining  in  varied  environments 
remain  the  major  uncertainty  in  application:;. 


5.1.1  Swell 

Ocean  swell  consists  of  long,  not  very  steep  waves  that 
have  run  faster  and  farther  (i.e.,  Cp^  = VglT  «VT”and  less  dissi- 
pation) than  their  shorter,  steeper  siblings  from  a storm  center. 

At  a distant  site,  typical  sv.ell  then  exceeds  sea  waves  in  the 
important  low-frequency  regime  when  local  winds  are  below  gale 
strength  (speeds  ^ 25  knots,  cf.  Fig.  2).  More  detailed  comparison 
of  swell  and  sea  will  be  made  in  Sec.  5.1.3. 

The  swell  has  a relatively  narrow-band  spectrum,  periods 
T~  10-20  s (Fig.  2).  Swell  wave-steepness,  6 = 2av/A,  is  far  be- 
low the  ideal  limit  (1/7),  or  even  the  observed  wind-wave  limits 
« 1/10  ).  A nominal  steepness  might  be  A ^ 1%,  in  which  case 
av  ~ 6 • T2  ^ 1-4  m (cf.  A ~ 160-640  m).  In  any  event,  from  Eqs. 

38  and  39a  we  have 


^surface  w a 10  ^ ^ Y * 


Swell: 


r T -4z/T  -z/(25  to  100) 

h , w a • T7  - e ^ e ' 

air  10 


(40) 


(41) 


with  T ~ 10-20  s and  a ^ 1 m. 

Thus,  for  example,  a magnetometer  of  sensitivity  10-1+  y 
could  register  swell-induced  0M  up  to  altitudes  of  order 

z ~ 2.4T2  ~ 200  m to  1 km  (42) 

max 


when  a ~ 1 m and  T ~ 10-20  s (and  to  depths  of  similar  order). 

We  note  that  in  the  case  of  shallow-water,  long-period 
swell  (Eq.  39b  above),  Eqs.  40  and  41  are  multiplied  by  a factor 
2>/d/T  < 1 and  the  exponent  in  Eq.  41  is  divided  by  the  same  factor. 
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Thus,  the  OM  field  becomes  smaller  and  less  extensive  (-*■  0 as  D 0) 
for  the  shallower  shelf  waters. 


5.1.2  Sea 

The  rms  OM  field  induced  by  a random  sea-wave  field  is  de- 
termined by  ensemble-averaging  Eqs.  38  and  39a  appropriately. 
Since  h « av,  this  means 


where  S(u>)  is  the  power  spectrum  of  wave  amplitudes  aw,  and  the 
factor  -g-  arises  from  time  averaging.  Sea  spectra  are  relatively 
broadband  (Fig.  2)  and  vary  with  wind  speed  (W) , fetch,  and  dura- 
tion. As  an  upper  limit,  one  can  consider  a fully  arisen  sea 
(long  fetch  and  duration,  so  only  W remains  as  a parameter)  and 
adopt  one  of  the  simple  phenomenological  spectra  that  have  been 
proposed.  In  general,  these  spertra  fall  exponentially  with  de- 
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creasing  y,  competing  with  the  opposite  exponentiation  in  h (Eq. 
39).  The  essential  result  is  a weakened  altitude  fall-off  for 
hj^  as  co»vared  with  h.  Thus,  the  spectral  forms  discussed  by 
Pierson  and  Moskowitz  (Ref.  19)  and  the  resulting  asymptotic  be- 
havior of  hi^  are,  aside  from  algebraic  factors, 


S (to) 


a 


h a 
rms 


m 

7 


» 


according  to  the  integral  approximation  given  below.  Here 
m=n/(n  + 2)  is  always  < 1 and  therefore  h^g  exponentiates  less 
rapidly  than  the  components  h in  Eq.  39a. 

For  more  specific  results,  we  will  adopt  the  generally  ac- 
cepted spectrum  recommended  by  Pierson  and  Moskowitz, 


S(u>) 


0.8l 

r _ u\ 

-BfT 

exp 

2 

nr-  - s 


» 


where  3 = 0.74.  The  older  Neumann  spectrum  has  been  used  in 
earlier  detailed  calculations,  which  are  summarized  into  simple 
formulae  at  the  end  of  this  section. 

Evaluation  of  the  integral  for  h^  at  the  surface  (z  = 0) 
is  elementary,  while  that  at  sufficient  altitude  can  be  approxi- 
mated by  the  general  asymptotic  formula 


f dwg(oj)e 
0 


■ i 

2tt  2 


-f(*>0) 

g(<Ve 


in  which  ui0  is  the  solution  of  f*  (w0)  = 0.  In  this  way  we  derive 
the  useful  simple  estimators 


Sea: 


h (z  = 0)  < 
rms 


30 


rms 


>1 

a 


W, 


K 


30 


(az)3  e y 


(43) 


(44) 


Ref.  19.  W.  J.  Pierson,  Jr.,  and  L.  Moskowitz,  "A  proposed 
spectral  form  for  fully  developed  wind  seas  based  on  the  similarity 
theory  of  S.  A.  Kitaigorodskii,"  J . Geophys , Res . . Vol.  69,  1964, 
p.  5181. 
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where  a e (3/2)“  gW2  = (7.65/WK)2  per  m and  is  wind  speed 

in  knots  (K).  Notice  that  Eq.  44  usually  obtains  at  reasonable 
altitude-  since  z > 1/a  < 10  m for  subgale  winds,  becoming  only 
~40  m fur  a whole  gale  (WK  ~ 50).  One  may  interpolate  roughly  be- 
tween z * 0 (Eq.  43)  and  z > 1/a  (Eq.  44). 

For  sea,  therefore,  a magnetometer  of  sensitivity  10~u  y 
could  register  the  induced  field  to  altitudes  (or  depths)  given 
roughly  by 


z 

~ W 2 • 

1 + -1  £n 

f w 

WK 

max 

2 K 

2 

L 30  j 

for  winds  from  a gentle  breeze  up  to  hurricane  strengths.  The 
height  (Eq.  45)  exceeds  1 km  only  for  strong  gales  and  decreases 
rapidly  for  lesser  winds. 

In  comparison  with  the  Person-Moskowitz  sea,  the  Neumann 
spectrum  S(u>)  « ui“6  exp  (-2g2/W2w2)  gives  a more  extensive  field. 
By  the  methods  above,  we  can  derive  formulae  corresponding  to 
Eqs.  43,  44,  and  45,  viz.: 


h (z  = 0)  < 
rms 


25 


rms 


f w 

f 11 

K 

TV 

A 

N 

< 

l 25  j 

(a'z)*  e-(“z>2 


where  a = 4g/W2  = (12.2/Wj^)2  per  n;  and  z 


max 


WK2  • [1  + *n(WK/25)]. 


5.1.3  Swell  versus  Sea 

Comparison  of  Eqs.  40,  41,  and  42  for  "typical”  swell  and 
Eqs.  43,  44,  and  45  for  a "full"  sea  shows  that  it  requires  a gale 
for  the  OM  field  induced  by  sea  to  exceed  in  size  and  extent  that 
induced  by  swell,  if  present.  That  is,  the  relative  steepness  (6 
up  to  10%)  of  sea  waves  is  more  than  offset  by  the  low-frequency 
content  of  ocean  swell,  unless  a gale  develops.  Notice  further 
that,  by  its  nature,  swell  is  widespread  in  contract  to  localized 
high-sea  conditions.  Thus,  swell  also  tends  to  dominate  sea  in 
horizontal  extent  over  the  ocean. 
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Figure  3 illustrates  the  swell-  and  sea-induced  OM  field 
magnitudes  given  by  Eqs.  41  and  44  as  functions  of  altitude. 
Curves  at  depth  display  similar  behavior,  except  that  h is  in- 
creased by  a factor  ~2k|z(  = 8|zj/T2  at  depths  |zj  ^ 1/k  « T2/4, 
where,  however,  exponential  decrease  also  becomes  dominant.  More 
detailed  curves  will  be  found,  for  example,  in  Ref.  4.  However, 
Fig.  3 is  sufficient  to  show  the  sizes  and  relative  regimes  of 
swell  and  sea  OM.  As  indicated  in  the  figure  caption,  Neumann 
sea  spectra  considerably  overestimate  the  OM  field  at  altitude 
relative  to  Pierson-Moskowitz  sea  spectra  estimates. 


5.1.4  Large  Wave  Crests 

We  have  so  far  C'  -idered  linearized  wave  theory.  Waves  of 
finite  height  (Ref.  20,  .ic.  250)  contain  higher  harmonics,  but 
these  are  of  relative  size  tea  = tt6  <<  1 even  for  the  steepest  sea 
waves.  Thus,  the  finit p amplitude  harmonics  introduce  negligible 
corrections  to  the  preceding  formulae.  However,  nonlinear  waves 
do  not  superpose  additively,  and  hence,  in  a high  sea,  chance  com- 
bination of  waves  can  form  individual  wave  crests  of  exceptional 
height. 

An  isolated  high  crest  will  produce  an  excess  contribu- 
tion to  the  OM  field  through  its  associated  mean  flow.  Naturally, 
this  cannot  be  treated  by  the  plane  wave  theory  as  above.  Instead, 
the  crest  OM  field  may  be  estimated  by  the  localized  flow  theory 
of  Sec.  4.2.  From  Ref.  20,  the  crest  possesses  a mean  flow  uQ  = 

tc2a2c  , e2<z  (z  £ 0,  c , = w/tc)  with  scale- length  L~  1/k,  so  that 
ph  ph 

its  magnetic  Reynolds  number  is  R*j  ~ paa2u  over  a depth  ~1/2k  and 
a horizontal  area  ~(ir/<)Lcr  where  Lcr  is  the  transverse  crest 
length;  i.e.,  Vol  ~ 7rLcr/2K2  in  Eqs.  30  and  31.  According  to  Eq. 
24,  this  flow  in  the  earth  field  induces  a local  OM  field  that  is 
~<a  <<  1 times  the  linear  wave  result  of  Eq.  38.  However,  at  dis- 
tances R » Lcr,  the  field  falls  off  as  1/R3  according  to  Eqs.  30, 
instead  of  exponentially  as  in  Eq.  39.  Further,  since  it  varies 
«a2  instead  of  linearly  with  amplitude  as  in  Eqs.  38  and  39,  it 
is  possible  for  an  exceptionally  high  crest  to  exceed  the  overall 
sea  in  the  higher  altitude  0M  field. 

The  crest  flow  induces  currents  whose  effective  dipole  mo- 
ment is  meff  ~ Heyoa2a)Lcr/16K2  according  to  Eq.  30b  with  the  above 
% and  Vol  expressions.  For  a deep  sea  (w2  « g<)  then  meff  ~ 
fycgHe  • (a2Lcr/4m<),  cf.  Eq.  38.  Hence,  at  ranges  R » Lcr,  we 
have  from  Eq.  30a 
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Fig.  3 OM  Field  Estimates  for  Pioi son-Moskowitz  Seas  and  Deep-Ocean 
Swells.  Dashed-line  (whole-gale)  sea  is  rarely  fully  arisen.  Note: 
Neumann-sea  OM  would  be  much  larger,  e.g.  the  fresh-breeze  OM 
would  approach  the  moderate-gale  curve  shown,  while  fresh-gale 
(37-kt)  OM  would  lie  near  the  whole-gale  curve. 
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Crest:  h 


1 a 
10  v 


aLcr 

4^F 


(46) 


where  the  factor  (aLcr/4KR3)  « 1/R3  now  replaces  the  exponential 
factor  in  Eq.  39. 

Again,  for  numerical  illustr  ition  we  assume  sensitivity 
10-4  y and  find  the  maximum  range 

R ~ 4(a2L  )3  T < 1 km  , (47) 

max  cr 


where  the  upper  limit  assumes  a long  T ~ 20  s,  an  Lcr  ~ 1/k  ~ 100m, 
and  a crest  height  a ~ 3 m which  might  occur  occasionally  in  a sea 
state  of  about  6 (Wj(  ~ 30  knots).  Thus,  an  isolated  high  crest  can 
sometimes  produce  an  OM  field  exceeding  that  of  swell  at  altitudes 
~ 1 km.  Notice  that,  since  the  crest  moves  at  speed  cph,  its  OM 
field  is  like  that  of  the  effective  dipole  moving  at  speed 
m/s  ~ T knots. 


5.1.5  Long  Waves 

We  turn  briefly  to  low-frequency,  long-wavelength  SW,  for 
which  the  nonstatic  limit  of  Eq.  34b  can  become  applicable  (cf. 

Fig.  1).  Aside  from  the  tides,  iong  ocean  waves  are  rare,  i.e., 
large  seismic  sea  waves  (tsunamis)  occur  approximately  once  a year. 
Nevertheless,  the  HM  effects  of  both  tides  and  tsunamis  have  been 
analyzed  and  compared  with  observation  by  Larsen  (Refs.  9 and  21). 
As  the  numerical  results  of  his  detailed  analyses  are  readily  esti- 
mated as  to  order  of  magnitude  from  our  theory,  we  shall  do  so  for 
completeness.  Our  estimates  will  not  include  the  finer  effects 
due  to  sediment  and/or  mantle  conductivities  (of  significance  only 
for  the  longest  waves),  for  which  see  Ref.  9 (Table  1).  Note  that 
his  "induction  parameter"  Q is  essentially  our 

A deep-ocean  tsunami  falls  under  the  nonstatic  regime  of 
Fig.  1.  Hence,  Eq.  34b  and  the  following  discussion,  together 
with  L~  1/k,  'uq  ~ agK/w  (see  Eq.  39b),  and  u>/k  5 cp^  <»  /gD  (Eq. 
35a) , gives  at  once 


Ref.  21.  J.  C.  Larsen,  "Electric  and  magnetic  fields  induced 
by  deep  sea  tides,"  Geophys.  J.  Roy.  Astron.  Soc. , Vol.  16,  1968, 
p.  47.  ~ " 
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(48) 


which  is  simply  the  tsunami-amplitude/ocean-depth  ratio.  For 
a ~ 0.2  m,  D ~ 5 km,  and  He~^  • 10s  y,  this  yields  0M  field  magni- 
tudes of  h ~ 2 y,  as  were  obtained  in  Ref.  9. 

A shallow-water  tsunami,  on  the  other  hand,  falls  into  the 
quasi-static  regime  of  Fig.  1.  Then  Eq.  39b  applies  and 


h ~ — avf)  y 


(49) 


Thus,  for  a *■“  0.2  m and  D ~ 200  m,  we  get  h ~ 0.6  y,  as  in  Ref.  9. 

The  tides  considered  as  free  waves  (e.g. , in  the  Atlantic 
or  Pacific  basins)  of  frequency  v < 2 • 10-5  Hz  have,  from  Eq.  35b, 


_h_ 

H 


KUr 


J^O 
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ph 
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D 


(50) 


the  same  as  for  the  deep-water  tsunami.  Modifications  due  to 
boundary  effects,  such  as  proximity  to  an  ocean  island,  are  cal- 
culated in  Ref.  21. 

All  of  these  OM  fields  are  very  large  scale,  comparable  to 
that  of  the  long  SW  (>100  km,  cf.  Fig.  1),  both  horizontally  and 
vertically  from  sea  bottom  to  the  air  above  the  sea  surface.  Thus, 
locally  they  appear  as  periodic  geomagnetic  noise  components. 


5.1.6  Bound  Waves 

It  is  possible  that  a high-frequency  system  of  bound  waves 
may  be  carried  along  by  a single  "carrier"  wave  (of  frequency  uic, 
speed  cc)  in  the  case  of  nonlinear  deep-water  waves  (cf.  Ref.  22). 
In  such  an  event,  the  bound  waves  have  a constant  speed  cc  and 
their  <a>(jc)  dispersion  relation  bacomes 


u>  = c k « t 
c A 


(uj  > w ) 
c 


Ref.  22.  H.  C.  Yuen  and  B.  M.  Lake,  "Reply  to  comments  of  G. 
Roskes,"  Phys . Fluids , Vol.  19,  1976,  p.  767. 
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which  is  steeper  than  the  dispersion  curve  for  the  usually  assumed 
free-wave  spectrum  ("sea"  in  Fig.  1).  Thus,  the  quasi-static  esti- 
or  of  Eq.  34a  remains  applicable,  leading  as  above  to  Eqs.  38 
and  39  times  the  factor  cc/c.  For  example,  for  deep  water  the 
bound-wave  spectrum  would  contribute  the  high-frequency  hydromag- 
netic  amplitudes 


h . » 

air 


x Eq.  39a 


As  a result,  for  a sea  spectrum  containing  bound  waves,  the  rms 
hydromagnetism  estimates  of  Subsec.  5.1.2  may  increase  as  follows: 
Eq.  43  for  the  surface  field  by  ft*  cc/W  and  Eq.  44  at  z > 1/a  by 
~ cc/W  (az/2)*/6,  provided  these  correction  factors  exceed  unity; 
otherwise  the  change  is  negligible.  However,  typical  sea  spectra 
have  little  contribution  at  wave  speeds  ^ wind  speed  (i.e.,  u> 
g/c  ^ g/W  in  Fig.  2),  so  that  one  expects  cc  < W.  Therefore,  the 
existence  of  bound  waves  would  have  little  effect  on  the  sea-in- 
duced hydromagnetism  curves  of  Fig.  3. 


5.2  INTERNAL  WAVES 

Buoyancy  waves  on  a density  stratification  (pycnocline) 
within  the  ocean  generate  OM  fields  that  are  usually  smaller  than 
the  SW-induced  fields  above.  However,  for  observations  with  a 
magnetometer  towed  in  the  vicinity  of  the  pycnocline,  the  IW-in- 
duced  magnetism  can  become  significant.  IW  magnetic  fields  have 
been  calculated  by  Beal  and  Weaver  (Ref.  7)  for  a sharp  pycnocline 
and  in  a recent  extension  of  their  work  (Podney,  Ref.  8)  for  the 
other  extreme  of  an  exponentially  stratified  ocean.  In  general, 
the  OM  field  due  to  IW  can  be  estimated  readily  much  as  for  S’’ 
above.  Thus,  the  quasi-static  Eq.  34a  (with  parameters  from 
Table  1)  gives 

h ~ Heuo  ^ ~ j ~ (in  y)  , (51) 


applicable  to  any  IW  of  supratidal  frequency  (cf.  Fig.  1). 

For  the  sharp  pycnocline  model  of  Beal  and  Weaver,  we  have 

u/k  ~ Hu/k  = He  , , where  c , = Yg(5p/p)//K  + k coth  <d  for  deep 
ph  pn 

ocean  (D  > 1/k)  from  Eq.  35b.  Boundary  conditions  (Eq.  1Gb)  to- 
gether with  directional  dependence  (cf.  Eqs.  39  and  following)  in- 
troduce a factor  (±i) , so  that  Eq.  51  yields  roughly 
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h ~ — ■A/  g 

pycnocline  10  W6^p 


ir  for  A <£  2ird  , (52a) 


for  A £ 2rd  . (52b) 


(Away  from  the  pycnocline,  h decreases  essentially  exponentially 
with  a scale  length  of  1/k.)  Numerically,  for  pycnocline  depth 
d = 100  m and  density  step  <5p/p  < 10“ 2 , one  finds  near  the  pycno- 
cline 


— ^ few-tenths  y/m 
a 


for  IW  of  a wavelength  A > 1 km,  and  K /A  less  for  A < 1 km  (cf. 
Beal  and  Weaver's  numerical  example). 

For  the  exponentially  stratified  model,  Podney  (Ref  8) 
chooses  to  express  h in  terns  of  the  IW  horizontal  flow  at  the 
surface  his  v,  (n=l,2,...).  Since  "u  ~ (<D/nm)  • v,  , Eq.  51 
yields  k0’”  Rtl’n 


D VkJ,l  ~ 1 1 


for  the  numerical  example  D = 1000  m,  = 1 cm/s  of  Podney,  in 
rough  agreement  with  his  detailed  calculations. 

The  more  detailed  behavior  of  the  0M  field  can  be  seen  from 
the  model  calculations  of  Beal  and  Weaver,  and  Podney.  However, 
Eqs.  52  and  54  display  clearly  the  essential  factors  that  deter- 
mine the  size  of  the  IW  magnetic  field. 

In  a stratified  fluid  there  can  also  exist  internal  shock 
waves,  or  hydraulic  jumps  (Ref.  23),  analogous  to  the  well-known 
SW  phenomena  (e.g.,  bores  and  solitary  waves,  cf.  Lamb,  Ref.  20, 
pp.  280  and  423).  Typical  of  earlier  observation  in  the  Straits 
of  Gibraltar  and  in  Massachusetts  Bay  are  the  detailed  measure- 
ments of  Hunkins  and  Fliegel  in  Ref.  24.  A nominal  example  has  an 


Ref.  23.  C.  H.  Su,  "Hydraulic  jumps  in  an  incompressible 
stratified  fluid,"  J.  Fluid  Mech. , Vol.  73,  1976,  p.  33. 

Ref.  24.  K.  Hunkins  and  M.  Fliegel,  "Internal  Undular  Surges 
in  Seneca  Lake:  A Natural  Occurrence  of  Solitons,"  J.  Geophys. 
Res. , Vol.  78,  1973,  p.  539. 
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amplitude  of  ~10m,  speed  of  ~§m/s,  and  time  scale  of  ~ 10  min, 
which  leads  via  Eq.  25  to  the  estimate 


h 


jump 


few  gamma 


for  this  transient  hydromagnetic  field. 


5.3  ACOUSTIC  WAVES 

Sound  in  the  sea  generates  a hydromagnetic  field  by  virtue 
of  the  interaction  of  the  acoustic  seawater  motions  with  the  earth's 
magnetic  field.  This  acoustic  OM  field  is  carried  along  with  the 
sound  wave,  thereby  radiating  through  th2  ocean  as  a pseudowave 
induction  field.  We  shall  refer  to  the  field  as  sonomagnetism,  to 
distinguish  it  from  the  magnetosonic  radiations  of  compressible 
MHD.  The  difference  is  that  sonomagnetism  propagates  purely  me- 
chanically via  the  acoustic  forces,  whereas  magnetosonics  involves 
in  addition  the  magnetic  action  back  on  the  fluid  motions.  In  the 
ocean,  that  reaction  is  very  small  (cf.  the  discussion  following 
Eq.  13) , and  so  magnetosonic  modes  collapse  to  ordinary  sound  waves 
(Ref.  16).  (Note  however  that  for  sonic  damping,  rather  than 
propagation,  magnetic  forces  become  more  important  than  viscosity 
at  infrasonic  frequencies.  That  is,  the  Hartmann  number  given  in 
the  discussion  following  Eq.  13  with  L ~ 1/k  = cs/co  becomes  % ~ 1/v 
for  the  ocean  parameters  of  Table  1,  so  that  magnetic/viscous  damp- 
ing > 1 for  v ^ 1 Hz.  This  happens  to  correspond  roughly  to  the 
transition  frequency  of  Eq.  55  below.)  Thus,  the  sound  wave  here 
will  be  taken  as  prescribed  by  the  usual  acoustic  equations  (Ref. 

25)  for  an  assumed  homogeneous  ocean. 

Sonomagnetism  was  discussed  theoretically  by  Kontorovich 
(Ref.  26)  for  plane  sound  waves  in  an  infinite  ocean,  with  nrief 
consideration  of  reflection-transmission  at  the  air-sea  interface. 

We  can  recapitulate  his  results  from  our  theory  and  discuss  in  ad- 
dition the  effect  of  finite  ocean  depth.  The  latter  is  necessary 
particularly  for  infrasonic.  frequencies,  since  in  that  regiir.t  the 
wavelength  can  become  comparable  to  depth  (cf.  Fig.  1).  and  alsc 
the  ambient  sonomagnetism  will  become  largest. 


Ref.  25.  H.  Lamb,  The  Dynamical  Theory  of  Sound  (1925),  2nd 
edition,  Dover,  1960. 

Ref.  26.  V.  H.  Kontorovich,  "Mag'.etohydrodynamics  of  the 
Ocean,"  Problems  of  MHD  and  Plasma  Dynamics  II,  2nd  Conf.  on 
Theoret.  and  Appl,  MHD,  Riga,  USSB,  1960. 
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For  a plane  sound  wave  in  infinite  ocean,  the  sonomagnetic 
pseudowave  has  identical  frequency,  wavelength,  and  velocity  as 
the  sound.  But  its  amplitude  (and  phase)  depend  markedly  on 
whether  the  frequency  (v)  is  above  or  below  the  rate  of  magnetic 
diffusion,  i.e.,  closer  to  the  quasi-static  or  nonstatic  limits, 
Eqs.  34a  or  34b,  respectively.  The  transition  frequency  is  (see 
the  discussion  following  Eq.  35c) 


“0  pocs2 

v0  E 27  = "27"*'  2HZ 


(55) 


for  the  ocean  parameters  of  Table  1^  With  L ~ 1 /k,  w = kcs,  and 
by  continuity  (Eq.  12b)  Tig  ~ (oi/ic)  6p/p  = csB  6P,  Eqs.  34  yield 
the  OM  amplitude  estimate 


h 


B • 6P 


(v  » vfl) 


(v  « vQ) 


(56a) 


(56b) 


#9  ^ 

where  3 = l/pcs4  is  seawater  compressibility  and  6P  the  acoustic 
pressure  amplitude. 


An  exact  expression  for  the  OM  field  is  easily  derived 
from  the  basic  Eq.  23  for  infinite  plane  waves,  viz., 


h = [ H - (H  • k)  -=7  ] • B • 6P 
~ ~e  ~e  ~ k-2  1 


vQ  - iv 


which  displays  the  dependence  on  field  (fig)  and  propagation  (jc ) 
directions  and  the  general  frequency  behavior 


1 + ~ 


(^-)2j  exp  i |tan- 1 


For  example,  a total  field  magnetometer  would  measure  amplitude 

|h  • Hg/Hg)  * Hg  • B • 6P  sin2  6/A  + (v/v0)2  (cf.  Eqs.  56),  where 

6 is  the  angle  of  sound  propagation  with  respect  to  the  earth  field. 
In  general,  the  phase  difference  between  the  sonomagnetic  and 
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Fig.  4 Ambient  Sonomagnetic  Spectra  for  Typical  Ocean  Noises  (Ref.  27)  and  a 
Chesapeake  Bay  Fish  (Croakers,  Ref.  28). 


acoustic  waves,  tan_1(v/v0),  is  90°  in  the  regime  of  Eq.  56a  and 
0°  in  that  of  Eq.  56b.  For  the  sonomagnetic  field  amplitude,  we 
can  rewrite  Eqs.  56,  using  Eq.  55  and  the  parameters  of  Table  1, 
as 


h~  10(S"86>/2°  . (22  + v2)-*  y , 


(57) 


in  terms  of  the  acoustic  intensity  measure,  S dB  re  1 Pa,  for  which 
ST  - 10S/2°. 


According  to  Eq.  57,  ambient  sound  levels  in  the  ocean 
(Ref.  27)  produce  the  sonomagnetic  "noise"  amplitudes  shown  in 
Fig.  4.  These  are  order-of-magnitude  estimates  only  and  for  the 


Ref.  27.  G.  M.  Wenz,  "Acoustic  ambient  noise  in  the  ocean: 
spectra  and  sources,"  J.  Acoust.  Soc.  Am.,  Vol.  34,  1962,  p.  1936. 
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infinite-ocean  plane  wave  approximation.  Thus,  it  applies  for 
widespread,  distant  sources,  with  no  specific  account  of  sea-air 
and  seabed  reflections,  channeling,  etc.  (See  Sec.  6 regarding 
strong,  local  sources  at  shallow  submergences.)  For  such  condi- 
tions, Fig.  4 shows  that  ambient  sonomagnetism  is  detectable  with 
magnetometer  sensitivities  of  ~10~4  y/r/viz  only  in  the  far  infra- 
sonic  regime. 


We  may  note  that,  for  a noise  souce  at  great  distances  (R) 
such  that  its  radiation  approximates  a plane  wave  at  the  magnetom- 
eter, the  sonomagnetic  transmission  loss-factor  is  the  square  root 
of  the  sonic  factor  (Ref.  28),  i.e.,  h « 1/R  for  spherical  spread- 
ing, l/t^R  for  cylindrical,  etc.  This  obtains  because  h « sound 
amplitude,  not  intensity,  in  Eq.  56. 

We  briefly  consider  boundary  influences.  The  sea-air  and 
seabed  interfaces  are  included  in  the  dispersion  relation  of  Fig.  1, 
from  which  we  see  that  the  transition  frequency  vo  of  Eq.  55  dif- 
fers little  from  D * » up  to  D ~ 500  m.  (In  shelf  waters  of 
D ~ 100  m,  vo  increases  to  ~ 5 Hz.)  In  general,  the  sonomagnetic 
field  can  be  constructed  as  a particular  integral  hp  of  our  basic 
Eq.  23  (with  *»  0)  plus  a complementary  solution  jic  chosen  to 
satisfy  the  boundary  conditions.  Now,  using  the  acoustic  wave 
equation  (4  + «2/Cg  )jjq  = 0,  one  can  see  that  the  Ansatz  hp  = 

-uo/(io)po  + (i>2/c  2)  • curl (H  x u0)  satisfies  Eq.  23,  as  well  as 

Eq.  22b.  Therefore,  jic  is  that  transverse-wave  solution  of  the 
homogeneous  vector  Helmholtz  equation  ((A  - iwpo)hc  = 0,  cf.  Eq.  23 
with  the  right-hand  side  zero)  determined  through  the  boundary 
conditions  (Eqs.  6a  and  6b)  by  the  value  of  Jip  at  the  interfaces. 

For  example,  at  the  air-sea  surface  (z  = 0)  we  find  for  the 
acoustic  modes  of  Eq.  35c  that 


hp(z=0)  = ^ sin  2a  • 


(H  • z)k  + (H  • k)z 

f*/  »v6 


i + v/v. 


(58) 


where  a = tan-1  <D/(£  + ^)tt  contains  the  depth  (D)  and  mode  (£)  de- 
pendence. In  terms  of  the  acoustic  cutoff  frequencies  v£  e 
(2 e,  + l)cs/4D  (cf.  Fig.  1 and  Eq.  35c  with  tc  -*•  0),  we  have  § sin  2a  = 

(vc/v)/l  - (vc/v)^  which  nas  a broad  maximum  at  frequency  v = Jl  vc. 

Thus  for  v ^ vc  to  v ~ several  vc,  Eq.  58  and  also  the  h^  de- 
termined through  it  are  comparable  with  the  infinite  plane-wave 
expression  given  following  Eqs.  56.  Therefore,  in  that  infrasonic 


Ref.  28.  R.  J.  Urick,  Principles  of  Underwater  Sound,  McGraw- 
Hill,  NY,  1975. 
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regime  the  complete  sonomagnetic  field  is  still  fairly  estimated 
by  Eqs.  56  and  57  when  air-sea  boundary  effects  are  included. 

The  unique  feature  of  the  sonomagnetic  pseudowave  vis-h-vis 
the  acoustic  wave,  however,  is  that  it  extends  into  the  air  above 
the  sea  where  the  transmitted  sound  is  negligible.  As  stressed  in 

Subsec.  2.2,  a surface  field  of  the  form  °=  e^1"  as  in  Eq.  58 

(see  also  Eq.  11a)  leads  to  an  aerofield  (cf.  Eq.  11b) 

iwt-iK  m:-kz 

h j = h(z=0)  • e <Z  « e . (59a) 

»v3xr 


That  is,  the  sonomagnetic  aerofield  looks  like  a surface  wave 
(pseudo  "ground  wave")  with  altitude  penetration  measurea  by 


Scale  height  = — ~ — — - — 

6 k (2,  + v>)lT 


(59b) 


near  the  optimum  frequency  v ~ Jl  v£.  Thus,  on  deep  ocean  and  for 
v in  the  0.1  to  1 Hz  band,  the  sonomagnetic  field  has  a scale 
height  of  ~ ktns.  This  is  much  greater  than  for  other  ocean-wave  OM 
in  that  frequency  regime  (viz.,  ^ 10  times  the  scale  height  for 
long  (20  s)  swells,  cf.  Eq.  41)  simply  because  infrasound  wave- 
lengths are  relatively  long  (Fig.  1).  Thus,  although  in  the  sea 
sonomagnetism  is  smaller  than  l*-generated  fields,  in  the  air  it 
can  remain  detectable  to  altitudes  beyond  the  limit  for  such 
fields.  For  example,  microseismic  sound  (Fig.  4)  might  in  princi- 
ple be  studied  free  of  surface-wave  noise  through  the  high-alti- 
tude sonomagnetic  field. 


5.4  SEACURRENTS 

Ocean  circulations  interacting  with  the  earth's  magneto- 
static field  have  long  been  recognized  as  sources  of  considerable 
electromagnetic  perturbations,  over  and  above  those  effects  in- 
duced in  a static  ocean  by  temporal  variations  in  the  earth's 
field.  In  1968,  Bullard  and  Parker  (Ref.  5)  reviewed  both  types 
of  oceanic  induction,  though  with  emphasis  on  the  latter.  The 
time-variation  effects  may  be  regarded  as  the  ocean  environment 
contribution  to  geomagnetic  noise  and  will  not  be  considered  fur- 
ther here.  The  other,  sea-current  effects  are  discussed  for  a 
simple  model  of  major  ocean  streams  in  Subsec.  4 of  Ref.  5,  with 
further  references  therein.  Later  references  are  cited  in  this 
subsection. 
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^Sea-current  OM  fields  can  be  estimated  from  Eqs.  24  and  25, 
h «?  ^ uqL  (y)  , (50) 

given  the  current  speed  and  scale,  uQ  and  L.  Of  particular  inter- 
est to  observers  has  been  the  electric  field,  E « curl  h/a  - pu0 x He, 
by  Eq.  7c  in  the  present  approximation  (h  <<  He,  jig  uniform),  or 
sometimes  only  that  part  measured  by  electrodes  moving  with  the 
seawater,  the  "geomagnetic  electrokinetic"  (GEK)  field,  j^  wcurl  h/a. 
Generally,  one  can  estimate  these  electric  fields,  using  [curl  h| 
h/L  and  Eq.  24,  as 

E ~ pHfiu0  ~ j 10-4  Uq  (V/m)  , (61) 


which  involves  only  the  current  speed,  not  scale. 

In  the  open  ocean,  Ekman  currents  driven  by  winds  of  ~ 10 
knots  may  have  surface  speed  Uq  ~ 10-1  m/s  and  can  extend  to 
depths  up  to  ~102m,  in  which  case 

h<l  V , E ^ pV/cm  . (62) 


Detailed  analysis  in  Ref.  29  confirms  this  OM  estimate.  Open 
ocean  tidal  currents  have  comparable  speed  and  therefore  electric 
field.  For  example,  Harvey  (Ref.  30)  measured  E as  ~ 0.014  pV/cm 
at  the  lunar  semidiurnal  tidal  frequency. 

Major  open  ocean  streams  can  have  higher  speeds  and  larger 
0M  effects.  For  example,  the  equatorial  countercurrent  (eastward) 
has  uq  up  to  1 m/s  and  depth  ~102m,  so  that 

h < 25  y , E < pV/cm  . (63) 


Ref.  29.  G.  A.  Burtsyev,  "Magnetohydrodynamics  of  sea  cur- 
rents," Fiz.  Atmosf.  Okeana,  Vol.  11,  1975,  p.  1084 

Ref.  30.  R.  R.  Harvey,  "Derivation  of  oceanic  water  motions 
from  measurement  of  the  vertical  electric  field,"  J,  Geophys.  Res., 
Vol.  79,  1974,  p.  4512. 
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(Cf.  Bullard-Parker's  simple  model  calculation  in  Ref.  5.)  Thus, 
Rommel  and  McCleave  (Ref.  31)  quote  their  unpublished  data  predict- 
ing up  to  0.46  pV/cm  in  the  Gulf  Stream  and  von  Arx  (corrected) 
figures  of  ~ 0.1 pV/cm  in  other  Atlantic  currents. 

Boundary  currents  on  eastern  margins  of  the  ocean,  e.g., 
the  California  current,  tend  to  be  sluggish  and  relative  shallow. 
Thus,  typically  uQ  ~ 1/10  m/s  and  L < 500  m,  so  that 

h < 10  y , E < pV/cm  (64) 


i 

I 

I 

I 

1 

i 

I 

I 

I 


I 

I 

I 


are  generous  upper  limits. 

Western  boundary  currents,  on  the  other  hand,  are  fast  and 
deep  enough  to  produce  considerably  larger  OM  fields.  For  example, 
the  Florida  current  of  the  Gulf  Stream  is  geostrophically  driven 
to  speed  uQ  of  ~ 1 to  3 m/s  with  a draft  of  <-  1 to  1.5  km,  so  that 

h £ 102  °r  3 Y , E ^ 1 pV/cm  . (65) 


All  of  the  above  sea  currents  are  very  large  scale,  with 
widths  typically  ~ 100  km.  Thus,  the  induced  OM  field  h is  corre- 
spondingly extensive,  comparable  to  the  scale  of  geologic  anoma- 
lies to  the  earth's  dipole  field. 

On  a smaller  scale,  one  may  consider  the  Langmuir  circula- 
tion. consisting  of  ? wind-induced  array  of  pairs  of  opnosing  heli- 
cal current  cells.  For  windspeeds  ^ 5 m/s,  the  cell  widths  are  5 
to  50  m (Ref.  32),  and  typical  current  speeds  are  ~ cm/s  (Ref.  33). 
Thus,  the  OM  field  forms  a periodic  spatial  pattern  of  wavelength 
(X)  twice  the  cell  width,  or 


h 10  2 t0_1  y with  A~10lto2m  . (66) 


Ref.  31.  S.  A.  Rommel,  Jr.,  and  J.  D.  McCleave,  "Oceanic 
electric  fields:  Perception  by  American  eels?"  Science,  Vol.  176, 
1972,  p.  1233. 

Ref.  32.  T.  Gammelsr«ld,  "Instability  of  Couette  flow  in  a 
rotating  fluid  and  origin  of  Langmuir  circulations,"  J.  Geophys. 
Res. , Vol.  80,  1975,  p.  5069. 

Ref.  33.  N.  E.  J.  Boston,  A.  Maratos,  J.  A.  Galt,  and  E.  D. 
Traganza,  "A  measurement  of  nearshore  Langmuir  circulation,"  Naval 
Postgraduate  School  Report  NPS-58Bb75041,  1975. 
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An  exact  calculation  (unpublished)  for  a simple  model  of  the  Lang- 
muir cells  finds  h reduced  by  a factor  ~ (cell-height/cell-width)2, 
which  can  be  an  order  of  magnitude  or  more.  Note  that  by  Eqs  11 
the  field  pattern  (Eq.  66)  extends  into  the  air  with 


Scale  height  - •—  < 15  m 


(67) 


so  the  Langmuir  circulation  OM  effects  are  correspondingly  con- 
fined  to  the  neighborhood  of  the  sea  surface. 


I 

! 


B 
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6.  ESTIMATES  OF  OM  INDUCTION 
BY  BODY-RELATED  SEA  MOTIONS 

The  size  and  extent  of  the  OM  fields  induced  by  seawater 
motions  related  to  the  passage  of  a body  through  the  ocean,  as 
listed  in  Eqs.  18e  through  18h,  are  now  estimated  from  the  theory 
of  Sec.  4.  As  indicated  there,  the  basic  Eq.  23  is  driven  by  a 
body  magnetization-translation  x V^)  interaction  and  bv  vari- 
ous geomagnetic- flow  (Hg  x jjo)  interactions.  Here  uq  wilr  be  sea- 
water disturbances  caused  by  the  body,  each  treated  separately 
(cf.  the  discussion  following  Eqs.  18)  and  in  isolation  from  the 
natural  sea  motions  considered  in  Sec.  5.  As  there,  clarity  and 
simplicity  are  chosen  over  precision,  but  some  detailed  analyses 
that  elaborate  our  estimates  will  be  cited. 

We  begin  (Subsec.  6.1)  with  the  magnetization-translation 
interaction,  simply  a mean-flow  effect,  assuming  a constant  body 
motion  and  magnetization  but  with  no  important  restrictions  on  the 
ocean  model  (uniform  conductivity  and  of  course  the  HM  approxima- 
tion (Eq.  4)).  Then  the  flow  effects  due  to  finite  body  displace- 
ment are  taken  up,  under  the  assumption  of  a homogeneous  and  in- 
compressible ocean.  We  first  consider  (Subsec.  6.2)  the  potential 
flow  around  a moving  body  at  large  Reynolds  number.  In  an  Oseen- 
like  approximation  for  the  OM  field,  the  potential  flow  effect 
will  be  seen  homologous  to  the  mean  flow  effect.  We  estimate  next 
(Subsec.  6.3)  wake-flcw  effects  in  uniform  or  stratified  ocean  via 
the  theorem  of  Subsec.  4.2.  Brief  consideration  is  given  (Subsec. 
6.4)  of  the  OM  effects  radiated  via  surface  wave  generation  by  a 
moving  body.  We  then  go  on  to  other  OM  pseudoradiations  that  arise 
by  virtue  of  the  inhomogeneity  an!  the  compressibility  of  the  ocean. 
Inhomogeneity  introduces  propagated  induction  via  internal  waves 
(Subsec.  6.5)  generated  by  a body  moving  through  stratified  ocean. 
Finally,  compressibility  (Subsec.  6.6)  allows  acoustic  waves  that 
induce  sonomagnetic  radiation  accompanying  sound  emissions  from  the 
body. 


6.1  MEAN  FLOW 

Consider  the  hydromagnetic  effect  that  arises  simply  from 
the  translation  through  the  ocean  of  the  magnetostatic  field  of  a 
moving  magnetized  vehicle.  That  is,  we  set  aside  until  later  any 
actual  seawater  disturbances  due  to  the  finite  displacement  of  the 
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body.  We  are  then  left  with  just  the  "disembodied"  magnetostatic 
field  Hjj  sweeping  through  the  seawater  and  thereby  inducing  mo- 
tional electrotaotive  forces,  eddy  currents,  and  magnetic  fields. 

Or  from  the  point  of  view  of  the  moving  body,  we  have  a mean  flow 
(-Vb)  of  seawater  through  the  magnetostatic  field  ) that  in- 
duces these  hydromagnetic  phenomena. 

The  mean  flow  effect  is  described  ^thematically  by  Eq.  23 
with  uQ  set  equal  to  zero.  We  assume  a constant  body  velocity  Vj, 
so  that,  in  the  steady  state,  3h/3t  ■ • 7)jj  on  the  left-hand 

s±de  of  Eq.  23.  The  spatial  dependence  in  Eq.  23  devolves  solely 
from  Hv  on  the  right-hand  side,  and  that  field  varies  algebraically 
with  distance  (R)  from  the  body  (r.f.  the  multipole  expansion  of 
magnetostatics).  Thus,  the  magnetic  Reynolds  number  for  the  mean 
flow  is,  by  Eq.  8,  R^  ~ paV^R,  and  the  first  (diffusion)  term  dom- 
inates on  the  left-hand  side  of  Eq.  23  for  ranges 


R « R 
v 


— V. 
pc  b 


(68) 


That  is,  the  quasi-static  approximation  (cf.  Eq.  24)  obtains,  and 
we  have  the  OM  estimate 


Si 


£ 


h ~ uoV.  RH.  ~ U°Vb  “b  (69) 

b 11  T? 

in  terms  of  the  magnetic  dipole  moment  of  the  body,  i.e.,  with 
Hb  ~ mjj/R3.  The  outstanding  feature  of  the  result  of  Eq.  69  is 
its  monopolar  range  behavior,  i.e.,  1/R2  falioff,  as  compared  with 
the  dipolar  1/R3  behavior  of  H^.  Of  course  this  is  restricted  to 
ranges  of  Eq.  68  and  to  an  infinite  ocean.  However,  we  see  below 
that  che  1/R2  behavior  persists  to  all  distances  in  the  wake  of 
the  body  and  in  a finite  ocean,  with  boundary  effects  mainly  intro- 
ducing complicated  angular  and  directional  dependence  (cf.  Eq.  71, 
below).  Thus,  Eq.  69  will  roughly  estimate  the  size  and  extent 
behind  the  body  of  this  nonisotropic  but  monopole-like  OM  field. 


Numerically,  Eq.  69  may  be  evaluated  by  using  m^  = /dt  ~ 
VM  • Mjj  where  is  the  volume  of  magnetic  matter  in  the  body  and 
its  magnetization.  (For  example,  a steel  shell  with  a thick- 
ness /diameter  ratio  6 would  have  ^ in  y - m3  ~ 106  • 6 • remanent 
magnetization  in  gauss  • body  displacement  4*r  tons.)  However,  it 
is  more  perspicuous  to  measure  h relative  to  Hb  from  Eqs.  68  and  69 
as 


*b 


D 

~ — for  R <<  R « 
R v 


400 


,Vb(in  kt)J 


km 


(70) 
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Thus,  h represents  a small  correction  to  that  increases  with 
distance.  It  appears  that  at  R~  Ry  (i.e.,  R^,  ~ 1)  we  might  have 
h ~ Hjj.  However,  there  the  second  term  on  the  right-hand  side  cf 
Eq.  23  becomes  important,  i.e.,  the  quasi-static  approximation 
breaks  down,  and  more  detailed  considerations  become  necessary. 


Tne  complete  analysis*  of  Eq.  23  for  a body  of  arbitrary 
magnetization  and  shape  includes  as  a particular  case  the  zero- 
displacement  mean  flow  effect  More  recently,  the  latter  case  for 
a point  dipole  in  a uniform  frow  has  also  been  giver,  by  Semenov 
(Ref.  34),  in  the  context  of  the  solar  wind  influence  on  the  earth 
dipole  field.  These  analyses  show  that  at  great  distances  (R  U-  Pv) 
the  magnetic  field  is  swept  back  to  form  a paraboloidal  wake  be- 
hind the  body,  analogously  to  Oseen's  solution  for  the  hydrodynamic, 
low  Re  flow  about  a sphere  (Lamb,  Ref.  20,  p.  613).  In  particular, 
our  analysis  shows  that  at  R > Rv  outside  the  wake  the  0M  field 
approaches  exponentially  the  negative  of  the  body  field,  i.e., 

£ -H^,  so  that  the  total  field  H = + h = 0 there.  On  the 

other  hand,  at  distances  R > Ry  within  the  wake  the  ON  field  re- 
tains the  1/R"-  behavior  of  Eq.  69,  plus  a 1/R  behavior  modified  by 
a small  factor  involving  the  downstream  angle  in  the  wake.  The  re- 
sult is  that  in  a downstream  "tail"  b now  remains  ~ 1/R^  as  in  Eq. 

69  at  all  ranges.  Thus  in  the  far  wake  (R  >>  Ry)  the  total  field 
is  essentially  just  the  OM  contribution,  i.e.,  h >>  H^,  so  that 
H = Hb  + h w h.  The  two  paraboloidally  separated  regions  of  dif- 
ferent behavior  of  the  OM  field  are  as  depicted  in  Fig.  5 of  Sub- 
sec. 6.2  (with  from  the  unified  analysis  equalling  m^  here). 

An  infinite  ocean,  or  plasma  in  Semenov's  case,  is  assumed 
in  the  above.  Otherwise,  the  only  noteworthy  restriction  is  that 
the  motion  and  the  moment  m^  be  constant  (cf . Ref.  34  for  oscil- 
lating or  precessing  m^) • 

Air-sea  boundary  effects  have  been  calculated  in  the  unified 
analysis  quoted  in  the  next  subsection.*  For  example,  for  a body 
with  vertical  moment  m^  moving  with  horizontal  velocity  at  depth 
d below  the  surface,  one  finds  the  results  given  in  Eqs.  79  and  80, 
below,  except  with  m0  there  replaced  by  m^.  That  is,  the  OM  field 
both  in  the  sea  and  in  the  air  above  has  the  form  (for  R <<  Ry) 


■w 


R2 


(71) 


*This  will  be  the  subject  of  a report  now  in  preparatior  by  J.  F. 
Bird. 

Ref,  34.  V.  S.  Semenov,  "Three-dimensional  dipole  in  a uni- 
form flow  of  conducting  fluid,"  Geomagn.  Aeron. , Vol.  15,  1975, 
p.  419  (APL/JHU  translation  2776). 
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in  which  {•••}  is  a complicated  vector 
unity.  Thus,  the  configuration  of  the 
{“•},  as  depicted  in  Figs.  6 and  7 in 
size  and  extent  of  the  field,  however, 
estimated  in  Eq.  69. 


function  of  magnitude  ~ 
OM  field  is  determined  by 
the  next  subsection.  The 
remains  roughly  (~ij)  as 


6.2  POTENTIAL  FLOW 

Now  consider  the  hydromagnetic  interaction  between  the  flow 
around  the  hull  of  a moving  body  and  the  local  magnetic  field.  In 
the  ocean,  the  hull  flow  may  be  taken  as  a potential  flow,  since 
the  Reynolds  number  (Re)  is  large  for  dimensions  and  speeds  of  in- 
terest here.  That  is,  Re  = Lu  p/n  « % * p/npo  by  Eq.  8,  which 
in  seawater  (Table  1)  gives  Re  ~ 10 11  • R^,  so  that  Re  >>  1 for 
any  % that  will  result  in  significant  OM  fields  (cf.  Eq.  24). 

The  local  magnetism  in  the  ocean  may  be  taken  as  a constant  uni- 
form earth  field  (cf,  the  discussion  following  Eq.  17),  ignoring 
any  body  magnetic  field.  That  is,  for  typical  body  magnetization 
one  has  Hj,/He  ~ nqj/HeR3  ~ (Mjj/He)  (Vfy/R3)  (cf.  the  discussion  pre- 
ceding Eq.  70),  which  is  small  everywhere  excepting  possibly  in 
the  immediate  vicinity  of  a thick-hulled  ferromagnetic  body.  But 
in  the  latter  case,  one  will  also  have  >>  h near  the  hull  so 
that  the  OM  field  is  swamped  by  the  body  field.  In  any  event,  we 
will  confine  our  attention  here  to  distances  outside  a few  body 
diameters. 

The  potential  flow  interaction  •• 'th  the  geomagnetic  field 
is  described  by  Eq.  23  with  the  right-  ia  d side  forcing  term 

curl(H0  x uQ)  -(H  • V)V  $ (72) 


where  <!>  is  the  flow  potential.  Now,  notice  that  the  body  term  in 
Eq.  23  treated  earlier  can  be  written 

-curios  x vfe)  = -(Vfe  • V)V  Ym  , (73) 

in  terms  of  the  scalar  magnetic  potential  to  the  body  field  H^. 
One  sees  at  once  that  the  potential  flow  effect  is  homologous  to 
the  mean  flow  effect.  In  general,  one  can  show  (Bird,  in  prepara- 
tion) that  the  potential  flow  OM  is  equivalent,  within  an  Oseen- 
like  approximation,  to  a "mean  flow"  OM  due  to  "translation"  Vj,  = 
& of  a "magnetization"  source  div  M = p(r),  when  the  potential 
flow  is  represented  by  a distribution  p(r)  of  sources  and  sinks. 
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In  particular,  a Rankine  ovoid  model  for  the  flow  gives  a 
very  simple  homology  relation  between  potential  and  mean  flow  ef- 
fects. The  ovoid  potential  at  sufficient  distances  R is  that  of  a 
source-sink  dipole  which  we  may  compare  with  a magnetic  dipole  po- 
tential: 


$ « s2a*  (V, 


(3>  • Z)  R 


(74) 


where  a and  s are  functions  of  body  length  L and  diameter  D 
(Ref.  14,  with  ict  m = our  s)  such  that  s2a3«»LD2/16  for  a slender 
body  (L  >>  D)  or  for  a sphere  (L  = D)  and  an  overestimate  by  at  most 
13%  for  intermediate  L/D  ratios.  From  Eq.  74  and  a commutation 
of  operators,  one  sees  Eq.  72  is  equivalent  to  Eq.  73,  only  with 
s2a3JJe  replacing  m^.  Therefore,  for  the  ovoid  and  at  R >>  a,  the 
potential  flow  OM  is  identical  with  that  arising  from  mean  flow 
Vb  past  a "hydromagnetic  moment" 


•Shm 


s2a3H 


LD2  H 
~e 


(75) 


Notice  that  depends  only  on  body  displacement  and  earth  field. 
The  purely  hydromagnetic  effect  (Eq.  75)  should  not  be  confused 
with  similar  excluded-volume  effects  of  magnetostatics. 

We  may  compare  mj,m  with  mj,  as  estimated  preceding  Eq.  70  to 
find  m^jj,  £ (Hg/%  • M^m^,  where  fM  is  the  fraction  of  total  body 
volume  that  is  magnetic.  Thus,  as  somewhat  of  a numerical  coinci- 
dence, a typical  ocean  vehicle  with  fM  < 10“ 3 and  M^,  10 3 gauss 

has 


m^  ^ m^  =>  potential  flow  OM  ^ mean  flow  OM  . (76) 


In  the  limit  of  a nonmagnetic  body,  of  course,  the  latter  is  ab- 
sent, but  the  potential  flow  OM  always  persists. 

The  combined  potential  and  mean  flow  effects  for  a moving 
magnetized  body  in  the  above  ovoid-dipole  approximation  can  now  be 
seen  to  be  equivalent  to  the  magnetization-translation  interaction 
for  an  effective  dipole 

”0  5 + Shm  * (77) 

Hence,  the  combined  OM  field  is  given  as  in  Eq.  69  with  Eq.  68, 
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y0V.ino 

h ~ — g2 — for  R « Rv  . (78) 

Again,  for  R £ Rv  the  field  h is  formed  into  two  paraboloidally 

separated  regions  as  depicted  in  Fig.  5,  with  the  OM  behavior  as 

described  in  Subsec.  6,1,  only  with  generalized  to  |q.  Thus, 

in  the  outer  region  h_  -*•  -Hq  exponentially,  where  Hq  is  the  field 

due  to  the  effective  dipole  j&q  of  Eq.  77,  so  that  the  total  field 

H = H + H,  + h = H - 1L  there  with  the  dipolar  1L  due  to  nv  . 

~ ~e  ~b  ~ ~e  ~hm  ~hm  ~hm 

(Note  that  the  body  field  is  still  cancelled  in  the  outer  region 

as  in  Subsec.  6.1.)  On  the  other  hand,  in  the  paraboloidal  wake, 

h^  retains  the  monopole-like  behavior  of  the  inner  region  given  in 

Eq.  78. 

Complete  details  of  the  OM  field  will  be  given  in  the  uni- 
fied analysis  to  be  published  (Bird,  in  preparation).  The  special 
cases  of  spherical  bodies  with  iJe  = 0,  ^ jj  (-JVj,  and  with 
Hgjjmb!!  (~Vjj)  were  treated  earlier  by  Ludford  and  Murray  (Ref.  35) 
and  by  Bois  (Ref.  36),  respectively.  These  are  included  as  degen- 
erate forms  of  our  analysis,  which  considers  nonparallel  iJe*  ^b* 
and  jgjj,  nonspherical  body,  nondipolar  and  bounded  media  (see 
below).  The  solution  was  obtained  by  an  Oseen-like  approximation, 
in  addition  to  the  usual  HM  and  incompressibility  assumptions  and 
constant  V^,  H^,  simplifications. 

The  air-sea  boundary  does  not  greatly  change  the  estimate 
of  Eq.  78,  but  mainly  the  detailed  configuration  of  the  OM  Held 
as  indicated  at  the  end  of  Subsec.  6.1.  To  begin  with,  the  effect 
of  the  air-sea  surface  on  the  potential  flow  is  accounted  for  by 
adding  to  the  ovoid  potential  in  Eq.  74  its  image  potential  (Ref. 
14).  Then  the  solution  of  Eqs.  23  and  22b  with  the  superposed  po- 
tentials is  found  that  satisfies  the  boundary  conditions  of  F-s.  6. 
Of  course  in  the  present,  finite-displacement  case,  one  must  also 
consider  the  body-sea  boundary.  The  latter  introduces  in  h terms 
~ 1/R4,  which  become  negligible  away  from  the  body.  If  the  body 
is  submerged  by  more  than  a few  times  its  dimensions,  this  body- 
sea  effect  is  unimportant  at  the  air-sea  boundary.  In  that  case, 
the  latter  independently  determines  the  far-field  effects  via  Eq. 
10b. 


Ref,  35.  G.  S.  S.  Ludford  and  J.  D.  Murray,  'On  the  flow  of 
a conducting  fluid  past  a magnetized  sphere,"  J.  Fluid  Mech.,  Vol. 
7,  1960,  p.  516. 

Ref.  36.  P.-A.  Bois,  "Influence  sur  quclques  types  de  champs 

magnfitiques  de  l'dcoulement  d'un  fluide  conducteur  autour  d’une 
sphfcre,"  J.  Mgc.,  Vol.  9,  No.  35,  1970. 
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Fig.  5 Hydromagnetic  Field  Regions  and  Wake  for  Body  Velocity  (Vb) 

in  a Medium  of  Magnetic  Diffusivity  (Miio)  (cf.  Table  1),  Leading  to 
Length  Scale  Rv  (Eq.  68).  The  OM  field  (h)  for  effective  magnetic  dipole 
(mo,  cf.  Eq.  77)  is  indicated  as  a function  of  range  (R). 


The  complete  solution  within  the  inner  region  and  paraboloi- 
dal wake  of  Fig.  5 (which  are  quite  extensive  in  the  ocean  — cf.  Ry 
in  Eq.  68)  is  detailed  for  such  a submerged  body  by  Bird  (in  prep- 
aration). For  example,  in  the  case  of  a vertical  effective  moment 
m0  and  horizontal  velocity  ^b,  one  finds  the  following  formulae: 

1,  Within  the  sea,  excepting  horizontal  ranges  comparable 
with  body  and/or  measurement  depths. 


h 


[R  cos  4> (1  + cos  0)  - 0 cos  <f>  sin  0 + $ sin  cf> ] 


poVbm0 

R2  5 (79) 


2.  In  the  air  overhead,  excepting  ranges  comparable  with 
body  depth, 


h 


R cos  4>  sin  0 - § cos  <f>  H-  <f>  sin  <)> 
1 + cos  0 


»*0Vo 

~~W~ 


. (80) 


% 

A 

I 


1 
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In  these  formulae,  R,  0,  4>  are  spherical  polar  coordinates  centered 
on  the  body,  with  the  polar  axis  vertically  tig,  azimuth  <|>  measured 
relative  to  direction  of  body  motion,  and  R,  6,  $ unit  coordinate 
vectors.  Field  configurations  calculated  from  Eqs.  79  and  80  are 
shown  in  Figs.  6 and  7. 

Aside  from  the  involved  angular  and  directional  dependence 
as  illustrated  for  this  example,  Eqs.  79  and  80  show  that  the  over- 
all size  and  extent  of  the  OM  field  due  to  combined  potential  flow 
and  mean  flow  effects  remain  roughly  as  estimated  by  Eq.  78. 


Fig.  6 Potential  Flow  and/or  Mean  Flow  Hydromagnetic  Configuration  in  the 
Sea.  Spatial  pattern  of  the  OM  field  parallel  to  motion  (h  ;|  Vb)  in  the 
horizontal  plane  centered  at  a body  (grid  units -y/paVb  mo/7). 
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Fig.  7 Potential  Flow  and/or  Mean  Flow  Hydromagnetic  Configuration  in  the 

Air.  Spatial  pattern  of  theOM  field  perpendicular  to  motion  (h  1 Vb)  in  a 
plane  at  altitude  15  grid  units  centered  above  a body  submerged  1 unit. 
(Grid  scale  1/2  that  of  Fig.  6). 


6.3  WAKE  FLOW 

Hydromagnetism  in  the  wake  of  a moving  body  will  be  consid- 
ered induced  by  a constant  uniform  earth  field,  neglecting  any  body 
field.  The  OM  effect  is  then  determined  by  the  character  of  the 
wake  flow,  i.e.,  whether  the  body  is  self-propelled  or  towed,  on 
the  surface  or  submerged,  in  a homogeneous  or  stratified  ocean, 
etc.  However,  estimates  of  the  OM  field  via  Subsecs.  4.1  and  4.2 
will  require  only  a general,  phenomenological  description  of  the 
wake-flow  magnitude  and  scale  (cf.  Eq.  25). 
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In  general  (cf.  Wu,  Ref.  37),  for  the  large  Re  ~ 1011  % 
indicated  earlier  (Subsec.  6.2),  the  wake  may  be  divided  into 
three  regions:  (a)  near  wake,  (b)  turbulent  far  wake,  and  (c)  lam- 
inar far  wake.  Farther  downstream,  a submerged  body  moving  in 
stratified  fluid  may  develop  a region  of  (d)  large  distant  eddies 
(Ref.  38).  In  addition  to  these  localized  disturbances,  of  course, 
there  can  be  radiated  disturbances  (IW,  SW,  AW),  that  will  be  con- 
sidered later  (Subsecs.  6.4  through  6.6). 

The  near  wake  extending  some  few  body  diameters  D down- 
stream, typically,  has  transverse  scale  ~D  and  flow  speed  ^ Vb  . 

The  local  OM  field  is  then  < HepaV^D  ~ by  Eqs.  24  and  25, 
while  at  distant  points  it  is  given  by  Eq.  31  as 

H yoV  D4 

h < (R  » D)  . (81) 

This  may  be  compared  with  the  potential  flow  effect,  Eq.  79  with 
m q ~ (l/10)D2LHe  and  L i*  D,  to  obtain 


near  wake  OM  I) 

potential  flow  OM  LR  R ’ 

which  is  small  since  R D. 

The  far  wake  develops  in  a complex  fashion,  particularly  in 
stratified  flow  (Refs.  39,  40,  and  41),  but  generally  the  wake- 
flow  magnitude  decreases  and  its  scale  increases,  eventually  col- 
lapsing vertically  in  stratified  fluid.  In  any  event,  the  turbu- 
lent far  wake  generates  random  OM  fields  of  approximately  zero 


Ref.  37.  T.  Y.  Wu,  "Cavity  and  wake  flows,"  Ann.  Rev.  Fluid 
Mech. , Vol.  4,  1972,  p.  243. 

Ref.  38.  H.-P.  Pao  and  T.  W.  Kao,  "Vortex  structure  in  the 

wake  of  a sphere,"  Phys.  Fluids,  Vol.  20,  1977,  p.  187. 

Ref.  39.  R.  J.  Hartman  and  H.  W.  Lewis,  "Wake  collapse  in  a 
stratified  fluid:  linear  treatment,"  J.  Fluid  Mech.,  Vol.  51, 
1972,  p.  613. 

Ref.  40.  R.  J.  Hartman,  "The  development  of  a partially 
mixed  region  in  a stratified  shear  flow,"  J.  Fluid  Mech.,  Vol.  71, 

1975,  p.  407. 

Ref.  41.  T.  W.  Kao,  "Principal  stage  of  wake  collapse  in  a 
stratified  fluid:  two-dimensional  theory,"  Phys.  Fluids,  Voi.  19, 

1976,  p.  1071. 
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mean,  but  of  rms  size  comparable  to  that  estimated  for  the  near 
wake.  The  laminar  far  wake  developed  by  relaminarization  further 
downstream  becomes  of  course  part  of  the  distant  potential  flow 
already  discussed  in  Subsec.  4.2. 

Finally,  consider  the  possibility  of  large  eddies  persisting 
farther  downstream,  as  indicated  by  the  experiment  of  Pao  and  Kao 
(Ref.  38).  in  stratified  fluid  of  small  Richardson  number 
(Ri  = N2D2/v^2  where  N is  buoyancy  frequency).  If  such  eddies 
are  produced  in  the  wake  collapse,  we  may  suppose  them  to  be 
slowly  rotating,  spheroidal  vortices,  with  vertical  size  ~D,  hori- 
zontal size  perhaps  an  order  of  magnitude  greater,  and  angular 
velocity  to  ~ N.  In  that  case,  Eqs.  24  and  25  yield  h < ND2,  or 
compared  to  the  near  wake  estimates  above, 


— .eddy,J)M  ^ M = fa 

near  wake  OH 


(83) 


which  is  small  for  Ri  <<  1.  Nevertheless,  since  these  putative 
eddies  are  very  distant  from  the  moving  body  and  large  compared  to 
it,  their  OM  field  configuration  may  retain  interest. 

For  example,  modeling  the  eddy  as  a spherical  vortex  of 
radius  a with  its  spin  vector  parallel  to  He,  one  can  use  the  de- 
tailed solution  given  in  Ref.  13,  Subsec.  9.9.  In  this  case,  the 
vector  jj  is  purely  horizontal,  azimuthally  directed,  with  magni- 
tude h = Hjpoco  a2 (sin  26/10)  • {(R/a)2  for  R €.  a,  (a/R)3  for  R ^ a} 
where  6 is  angle  from  vertical.  Thus,  the  size  of  h is  as  esti- 
mated above  (Eq.  83),  only  now  for  a sphere  and  with  a factor  s:l/10 
introduced  by  the  spherical  boundary  condition,  viz.,  at  most  h < 
wa2/40  at  R = a.  The  fall  off  « l/R3  outside  the  vortex  is  as  ex- 
pected from  the  theorem  of  Subsec.  4.2.  The  detailed  field  con- 
figurations in  horizontal  planes  passing  through  and  near  the 
spherical  vortex  are  illustrated  in  Fig.  8 (plane  at  a/5  above  or 
below  vortex  center)  and  in  Fig.  9 (plane  at  2a/5  above  or  below 
vortex  edge).  Modifications  due  to  air-sea  boundary  effects  are 
small,  as  long  as  the  eddy  is  submerged  more  than  a few  radii. 
However,  we  note  that  this  vertical-spin-vertical-field  example 
is  exceptional  in  that  its  OM  field  in  the  air  above  the  sea  is 
identically  zero  (cf.  following  Eq.  28). 

More  general  cases  of  arbitrary  eddy  shape,  spin,  and  field 
direction  will  have  field  configurations  outside  the  eddy  given  by 
the  theorem  of  Eq.  29. 
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Fig.  8 Vortex  Hydromagnetism  in  Plane  Through  Spherical  Eddy:  Spatial 
Configuration  of  Total  OM  Field  Over  Planes  of  Rotational  Motion 
Located  1/10  Eddy-Diameter  Above  or  Below  Eddy  Center. 


h 

t 


Fig.  9 Vortex  Hydromagnetism  Above  or  Below  Spherical  Eddy:  Spatial 
Configuration  of  Total  OM  Field  in  Planes  or  Rotation  at  1/5  Eddy- 
Diameter  Away  from  the  Eddy  Top  or  Bottom  Boundaries. 
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6.4  SURFACE  WAVES 

The  OM  field  due  to  body-generated  SW  is  estimated  using 
Subsec.  4.3,  somewhat  as  we  did  for  natural  SW  at  the  beginning 
of  Subsec.  5.1.  In  the  present  case,  however,  for  constant  body 
motion  V^,  a stationary  (in  body  frame)  wave  system  will  develop 
such  that  the  frequency  to  of  each  Fourier  component  with  wave  vec- 
tor k satisfies 


to  = k • V,  = kV,  cos  6 
~ ~b  b kV 


(84) 


Equation  84,  together  with  the  dispersion  relation  Eq.  35a,  delim- 
its the  permissible  waves,  as  detailed  below.  However,  for  any 
ordinary  body  speed  V},,  Eq.  84  alone  will  suffice  to  estimate  the 
OM  field.  That  is,  Eq.  84  plots  as  a line  of  slope  -45°  and  ordi- 
nate intercept  cos  6 Ky  in  Fig.  1,  which,  as  long  as  < 102m/s, 
will  always  intersect  the  SW  dispersion  curves  well  within  the 
quasi-static  regime.  Hence,  Eq.  34a  applies  with  u aw/tanh  kD 
and  L ~ (tanh  kD)/k  (Subsec.  5.1)  to  give  *h  ^ H^yoaw/ic.  Then  Eq. 

84  yields  at  once 

h ^ HepoV^  cos  0^  “a  , (85) 

as  an  upper  limit  on  the  OM  amplitude  in  terms  of  plane  SW  ampli- 
tude a. 


A limitation  for  Eq.  84  to  be  compatible  with  the  disper- 
sion Eq.  35a,  which  says  that  u < io^D,  is  the  requirement 

cos  e^v  s:  v^gD  • (86) 


For  most  usual  speeds  and  depths,  V^2/D  s g so  that  this  condition 
is  fulfilled  for  all  angles  0Ky.  Furthermore,  ordinarily  Vk2/D  <<  g 
so  that  w <<  lo'gD,  i.e,,  the  deep-water  approximation  is  applicable 
(cf.  Eq.  87,  below).  Only  for  very  fast  motion  in  very  shallow 
water  (V^2/D  > g)  is  the  above  restriction  important,  in  which 
case  it  limits  the  range  of  0Ky  (see  below) . 

In  the  ordinary  case  (Vfc2/D  <<  g) , more  accurate  formulae 
for  the  OM  field  analogous  to  those  in  Subsec.  5.1  are  readily  ob- 
tained. The  only  difference  here  is  that  Eq.  84  combined  with  the 
deep-water  dispersion  relation  u,2  w Kg  determines 
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< f* 


V,  * cos2  9 

D 


and 


ta 


tcV 


V,  cos  0 „ 
b kV 


(87) 


Notice  these  give  lower  bounds  (k  = g/V^2,  m - g/V^)  attained  for 
plane-wave  propagation  along  the  direction  of  body  motion  (6Ky  = 0) , 
The  first  element  of  Eq.  87  shows  that  the  body-generated  SW  are 
short  compared  to  depth  (<D^gD/V^2  >>  1)  and  the  deep-water  ap- 
proximation is  quite  good.  Also,  the  quasi-static  approximation 
of  Eq.  34a  (u/k2  <<  1/po)  is  verified  from  Eq.  87  to  be  adequate 
for  V]j  < (g/po)1'  sec  0Ky,  which  is  assured  in  the  ocean  (Table  1) 
for  Vf,  s;  102  m/s  (cf.  above).  Therefore,  just  as  in  Eq.  38  and 
the  discussion  following  Eq.  39,  only  with  '>  therein  given  by  Eq. 

87,  we  have 


surface 


ta 


1 *. 
16aVb 


{i, 


V 


1 tT 

l6aVb 


in  y 


(88) 


where  {I,  0K,  0y}  = cos(0K  - 8y)^sinz  I + cos2  I cos2  0K  in  terms 
of  dip  angle  I and  compass  angles  f*K,  0y  of  wave  propagation  and 
body  motion,  respectively,  so  } ©KV I = I etc  " evl*  Equation  88  gives 
the  OM  plane-wave  component  at  the  air-sea  interface  in  terms  of 
the  plane  SW  amplitude  “a.  The  altitude  (or  depth)  variation  is  ex- 
ponential (or  predominantly  so)  corresponding  to  Eq.  39  and  follow- 
ing. However,  in  the  present  case  one  sees  from  Eq.  87  that 

V 2 

1 v 

vertical  scale £ (89) 

Kg 

is  typically  small  (i.e.,  the  body-induced  SW  are  short)  so  that 
the  OM  field  is  essentially  restricted  to  the  near-surface  region. 
There  Eq.  88  will  give  a fair  estimate  of  the  OM  field. 

In  the  fast-shallow  situation  (V^2/D  > g) , we  saw  that  the 
permissible  wave  solutions  are  restricted  to  the  directions 
cos  0Ky  s:  v^gD/V^  (cf.  also  Ref.  20,  p.  440).  However,  Eq.  88  with 
this  restriction  remains  a fair  estimate,  as  does  Eq.  89  well 
within  the  allowed  directions.  For  0^  approaching  the  cut- 
off angle  0C,  the  scale  by  Eqs.  35a  and  84  becomes  1/k  = 

1//0kV  “ ec  -*■  “,  but  the  OM^amplitude  remains  bounded  in  accord 
with  Eqs.  85  and  86,  viz.,  h ~ Hepov^gD  sf  (cf.  Eq.  39b). 

For  a simple  illustration,  consider  the  one-dimensional 
(0K  ■=  0V  = 0)  SW  generated  by  a cylinder  of  radius  b moving  at  sub- 
merged depth  d in  deep  water.  According  to  Dagan  (Ref.  42),  linear 

Ref.  42.  G.  Dagan,  "Free-surf ace  gravity  flow  past  a sub- 
merged cylinder,"  J.  Fluid  Mech. , Vol.  49,  1971,  p.  179. 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL  MAR>LAND 


wave  theory  is  valid  ford  & 5b,  in  which  case  the  wave  amplitude 
is  'a  = 4irb2Ke  K<*  (Ref.  20,  p.  412).  Then  Eq.  88  with  Eq.  87  gives 


h « 


^b2g 


4V. 


exp 


y b2 


(d  £ 5b) 


(90) 


which  depends  on  cylinder  size  b and  depth  d as  one  might  expect. 
The  body-speed  dependence  is  such  that  the  0M  field  is  a maximum 
for  Vjj  = /2gd,  giving  h^  & b^/’/d  y < i|b3/2  y for  d ?>  5b,  from 
Eq.  90.  For  moderate  depths  (2.5b  ^ d < 5b),  the  same  formulae 
hold  with  d replaced  by  a larger  effective  depth;  for  shallower 
depths  (0  < d < 2.5b),  noulinearities  modify  the  SW  significantly 
(Ref.  42). 

Realistic  two-dimensional  wave  systems  require  Fourier  in- 
tegration of  Eq.  88  to  determine  the  OM  field  Ji.  Of  course,  h 
may  be  estimated  from  a knowledge  of  the  dominant  SW  amplitudes. 

The  intricacies  of  vave  patterns  (Ref.  43)  will  bo  reflected  in 
correspondingly  complex  configurations  of  the  OM  field.  However, 
perhaps  the  major  difference  from  one-dimensional  models  is  that 
the  SW  amplitudes  "a  fall  off  with  distance  from  the  source,  except 
possibly  in  any  singular  directions. 

For  example,  Kelvin  ship-wave  patterns  as  sketched  in  Fig. 

10  consist  of  transverse  ( /2/ 3 £ cos  £ 1)  and  lateral  (0  < 
cos  0kV  « /2/3)  systems  of  equal  amplitude  (Ref.  20,  p.  436).  How- 
ever, the  transverse  waves  being  longer,  they  have  larger  RM  ('“u/k 
« w/k  « cos  0Ky,  by  Eq.  87)  and  therefore  tend  to  contribute  more 
to  the  size  (and  extent  away  from  the  surface)  of  the  OM  field. 
Mathematically,  one  sees  in  Eq.  88  that  the  angle  factor  {I,  0K,  0v) 
is  larger  for  transverse  plane  waves  than  for  lateral,  excepting 
east-west  ship  headings  (0y~±n/2)  in  the  tropics  (I  ~ 0).  Partic- 
ularly in  polar  regions,  or  for  more  north-south  ship  headings 
elsewhere,  the  transverse  waves  dominate  in  the  0M  effect,  giving 


h 


sin  I 


Y 


(91) 


Of  course,  the  plane-wave  amplitudes  “a,  and  thus  h,  must  be  deter- 
mined from  the  usual  two-dimensional  circular  wave  solutions.  How- 
ever, well  within  the  wake  and  far  downstream,  the  transverse  sys- 
tem looks  like  plane  waves  of  slowly  decreasing  amplitude, 


Ref.  43.  C.  Hunter,  "On  the  calculation  of  wave  patterns," 
J.  Fluid  Mech.,  Vol.  53,  1972,  p.  637. 
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a = R 2 


«w  ~ °> 


(92) 


at  horizontal  ranges  R in  given  direction  (Ref.  20,  p.  436),  and 
the  OM  field  falls  off  correspondingly.  We  remark  that  in  the 
cusp  directions,  the  wave  height  decreases  more  slowly, 

a « Rc"2  (cos2  6kV  = | , 0kV  « ±35°)  (93) 

for  a (not  its  Fourier  amplitude  a)  as  a function  of  uistance  Rc 
along  a cusp  (Ref.  44).  Hence,  h can  fall  off  slowly.  However, 
more  elaborate  analysis  is  necessary  to  include  the  singular  direc- 
tions (Eq.  93) 


— 1 h~1//R 


Vb  9kV 


Fig.  10  Kelvin  Ship-Wave  Hydromagnetism:  h Contributed  Largely  by  Lateral 
Waves  (bold  wavefronts,  propagation  vector  k.)  From  Body  Surface 
Movement  (Vb),  with  Asymptotic  Range  (R)  Dependence  Indicated 
(cf.  Eq.  90).  Inset  shows  direction  angles  of  text  (cf.  Eqs.  84  and  88). 


Ref.  44.  F.  Ursell,  "On  Kelvin's  shipwave  pattern,"  J.  Fluid 
Mech. . Vol.  8,  1960,  p.  418. 
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6.5  INTERNAL  WAVES 

In  stratified  ocean,  a moving  body  generates  IW  whose  OM 
effects  can  be  estimated  as  for  the  SW  above.  IW  motion  is  in 
general  three  dimensional  and  vertically  inhomogeneous,  but  its 
horizontal  dependence  may  usefully  be  expanded  i-'  plane  waves  to,  jc. 
Then  for  constant  body  motion  V^,  again  one  has  a wave  system  sta- 
tionary in  the  body  frame  that  satisfies  Eq.  84  as  before.  Since 
the  IW  curves  in  Fig.  1 lie  considerably  below  the  SW  ones,  the 
compatibility  condition  between  Eq.  84  and  the  IW  dispersion  Eq. 

35b  will  be  more  restrictive  than  in  the  case  of  SW  (see  below). 

By  the  same  token,  however,  the  allowed  IW  a fortiori  fall  further 
than  did  the  SW  within  the  quasi-static  regime  of  Fig.  1.  Hence, 
Eq.  34a  with  Eq.  84  gives  as  above  immediately 

h(z)  < HeyoVb  cos  0kV  aj(z)  (94) 


for  the  OM  amplitude  due  to  a plane  IW,  in  terms  of  its  azimuth 
0Ky  in  the  horizontal  plane  corresponding  to  vetical  coordinate 
z and  its  amplitude  a^z). 

The  restriction  imposed  by  the  IW  dispersion  relation  can 
be  seen  from  Fig.  1,  again  imagining  there  the  -45°  line  (Eq.  84). 
It  intersects  the  IW  curves  only  if  it  does  not  exceed  the  limit 
for  long  IW,  for  which  the  thin  thermccline  of  Eq.  35b  is  most  apt, 
giving  a)  < K-/g(6p/p)d.  Therefore,  from  Eq.  84, 


cos  s /g(<5p/p)d  ~ few  m/s  (95) 

for  typical  pycnocline  depth  d and  density  change  5p/p  (Table  1) 
in  ocean  of  depth  D >>  d.  (For  finite  D,  Eq.  95  is  reduced  by 
vl  - d/D.)  Thus,  only  for  slowly  moving  bodies  or  for  very  deep, 
strong  pycnoclines  is  there  no  serious  restriction.  Otherwise, 

Eq.  95  delimits  permissible  IW  propagation  angles  0Ky  for  all 
Vfc  > /g(<5p/p)  (cf.  Ref.  20,  p.  416  for  the  one-dimensional  case, 
0KV  = 0) • 


To  illustrate  the  more  detailed  behavior,  we  may  consider 
various  .limiting  cases  that  arise  as  V^,  cos  0Kv  is  decreased  from 
its  limit  (Eq.  95).  One  sees  from  Fig.  1 that,  as  the  -45°  line 
of  Eq.  84  is  lowered,  it  intersects  the  IW  curves  initially  in  the 
long-wave  (1/k  ^ d)  region;  then,  for  cos  0Ky  somewhat  below 
the  limit  of  Eq.  95,  in  an  intermediate  wavelength  (e  <<  1/k  ^ d) 
region,  provided  the  pycnocline  structure  so  allows  (i.e.,  width 
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e <<  depth  d) ; and  ultimately  in  the  short-wave  (1/k  < e)  region 
where  the  IW  dispersio-i  curve  becomes  bounded  by  the  buoyancy  fre- 
quency, i.e.,  w2  £ N2  ~ g (<5p/p)/t  (cf.  Ref.  18). 


For  the  long  waves  with  cos  0Kv  w /g(<5p/p)d,  i.e.,  near 
the  cutoff  angles  (6C)  indicated  by  Eq.  95,  the  dispersion  Eq.  35b 
together  with  Eq.  84  gives  that 

k2  « — tan  6 • 1 6 , - 6 [ , (0  0 ) , (96a) 

Dd  c 1 kV  c'  kV  c 

again  for  D >>  d.  Thus,  k (and  n)  + 0 as  6KV  -*•  0C,  while  Eqs.  94 
and  95  with  the  values  of  Table  1 yield 


h(z)  /(dpTpTd  • a][(z)  Y , (6kV  + (96b) 


as  the  absolute  upper  bound  on  the  OM  air-ilitude. 

For  the  solutions  with  cos  GkV  somewhat  smaller,  but  not 
so  much  as  to  vitiate  the  thin  pycnocline  approximation,  say 

/g(6p/p)e  <<  V}j  cos  0Ky  < /g(6p/p)d/2  (assuming  e <<  d) , the  dis-  

persion  relation  is  the  Kd  ^ 1 form  of  Eq.  35b,  i.e.,  to  /xg(6p/p)/2. 
This  together  with  Eq.  84  determines  (cf.  Eq.  87) 


« and  M <»  ^,tL-eJpb  ■ , 

2Vb2  cos2  'vv  2Vb  cos  vv 


(97a) 


which  with  the  above  limits  on  cos  0^  span  the  wavelength  and 
frequency  ranges 


e <<  — ^ d and  SzTld  N ^ to  <<  N 


(97b) 


indicated  earlier.  Observe  from  Eqs.  97  that  the  quasi-static  ap- 
proximation of  Eq.  34a  (coL2  ^ oo/k2  <<  1/po)  is  assured  for  d << 

(2/p2o2g(6p/p)) 1^3,  which  is  at  least  ~ 104  m. 

This  range  of  cos  0Kv  leading  to  Eq.  97b  covers  the  bulk 
of  the  admissible  IW  for  a thin  pycnocline  (e  <<  d).  For  such 
waves.  Beal  and  Weaver's  work  (Ref.  7),  which  led  to  Eq.  52  above 
except  with  Cp^  = oj/k  = cos  0kV  here  by  Eq.  84,  gives  the  (±50%) 
estimate 

h ~ cos  0KV'aIw  jr  / (6p/p)/<  a^.  y (at  pycnocline)  (97c) 
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using  Eq.  97a.  Away  :rom  the  pycnocline,  h falls  essentially  ex- 
ponentially with 


f 

j 


vertical  scale 


(97d) 


In  Eqs.  97c  and  97d,  k varies  according  to  Eqs.  97a  and  97b. 

Eventually,  IW  with  cos  0Ky  /g(6p/p)e  enter  as  higher 
inodes  if  the  pycnocline  is  sharp  (e  <<  d) , or  in  place  of  the 
above  "thin-thennocline"  mode  if  it  is  diffuse  (e  ~ d).  In  either 
event,  these  waves  have  scale  length  s;  pycnocline  thickness,  so  de- 
tails of  the  stratification  assume  importance.  However,  one  can 
roughly  estimate  the  OM  field  from  the  general  properties  of  the 
IW  modes.  Their  spatial  variation  is  dominated  by  the  shorter  of 
the  two  lengths  given  in  Phillips  (Ref.  18,  Sec.  5.2),  the  hori- 
zontal  ~ 1/tc  for  to  ~ N and  the  vertical  fcz  ~ ih/^^/to2  - 1 ~ 
u/Nk  for  on  <<  N.  Thus,  one  can  write  for  the  scale  length  in  the 
mode  with  k = Kn, 


to 

^n  ~ nIT-  (mode  n)  . (98a) 

n 

Then  from  the  quasi-static  Eq.  34a  with  L ~ £n  there,  one  has 
directly 

h(z)  < Hepo£n  ux(z)  » ■£  Uj(z)  Y (98b) 


for  an  upper  limit  on  the  OM  amplitude  due  to  the  IW  mode  n with 
velocity  amplitude  uj(z). 

The  direction  angle  0Ky  for  these  higher  modes  by  Eqs.  84 
and  98a  satisfy 


cos  0 


tcV 


V.K 
b n 


~ N 


£ (mode  n) 


(99) 


I 

I 


This  represents  a quantization  of  0Ky  that  can  have  import  for  the 
extent  of  the  OM  field,  as  discussed  later.  For  the  OM  size,  how- 
ever, Eq.  99  simply  reproduces  the  estimate  of  Eq.  98b  with  uj  -*• 
to  aj  , using  Eq.  94  corrected  by  a factor  ~ to/N  to  include  the 
cases  &z  « encountered  here.  Note  that  by  assumption  above 

Vb  cos  0Ky  ^ /g(6p/p)e,  which  by  Eq.  99  implies  < e as  we  have 
indicated.  Also  notice  that  in  the  limiting  extreme  of  an  infinite 
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ocean,  in  which  case  d -*■  ®,  Eq.  95  no  longer  restricts  but  instead 
is  replaced  by  Vj,  cos  6<y  < Ne. 

To  determine  the  full  extent  and  configuration  of  the  OM 
field  requires  a detailed  analysis  equivalent  to  superimposing 
(via  multiple  Fourier  integration  and  with  proper  selection  to  fit 
the  boundary  conditions,  e.g.,  Eq.  10b)  the  various  IW  solutions 
we  have  estimated.  In  special  cases,  however,  some  estimate  of 
the  range  behavior  is  easily  obtained.  For  example,  in  a thin- 
thermocline  approximation  the  dominant  IW  displacement  at  the  in- 
terface behaves  asymptotically  like  a plane  wave  of  amplitude 


at  large  horizontal  ranges  p,  according  to  Hudimac  (Ref.  45,  Eq. 
44).  Then  the  OK  field  behaves  similarly  by  Eqs,  96b  or  97c.  At 
the  other  extreme,  in  an  infinitely  thick  constant-N  thermocline, 
the  IW  flow  velocity  behaves  asymptotically  like  a plane  wave  of 
amplitude 


at  large  radial  distances  R,  according  to  Miles  (Ref.  46,  Eqs.  6.12 
and  6.13).  In  this  case,  Eq.  98b  (with  £n  ~ Miles'  1/k  = Vb/N) 
gives  similar  range  behavior  for  the  OM  field.  More  generally, 
however,  the  IW  mode  structure  has  to  be  taken  into  account. 

The  quantization  of  the  IW  into  discrete  modes  can  lead  to 
singular  behavior  in  preferred  directions.  Consider  the  field  at 
large  distances  in  or  near  the  pycnocline,  which  corresponds  to 
small  k in  the  Fourier  inversion,  and  therefore  to  large  above. 
Thus,  the  higher  modes  for  which  >>  lz  contribute  largely  to 
this  asymptotic  field.  For  such  inodes  the  scale  length  is  £n  £z. 
Now  lz  is  determined  by  the  boundary  conditions  at  the  air-sea  sur- 
face and  seabed  (Ref.  18).  Therefore,  from  Eq.  99  the  mode  angle 
0Ky  is  also  essentially  fixed,  viz.. 


cos  9^,  * cos  0(n)  ~ -y-  &z(n)  (mode  n) 

b 


(ICO) 


Ref.  45.  A.  A.  Hudimac,  "Ship  waves  in  a stratified  ocean," 
J.  Fluid  Mech..  Vol.  11,  1961,  p.  229. 

Ref.  46.  J.  W.  Miles,  "Internal  waves  generated  by  a hori- 
zontally moving  source,"  Geophys.  Fluid  Dyn. , Vol.  2,  1971,  p.  63. 
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As  a consequence,  each  higher  mode  yields  via  Fourier  inversion  an 
asymptotic  field  contribution  that  is  predominantly  in  the  recipro- 
cal direction 

<Kn)=f-e(n)  , (101) 

where  <}>  is  the  field-point  horizontal  azimuth,  measured  from  the 
downstream  (-Vb)  direction.  (Briefly,  Eq.  101  arises  from  Eq.  100 
as  follows  : in  cylindrical  polar  coordinates  p,  4>,  z with  z axis 

vertical,  jj(p,<J>,z)  = f djce:i'Kp  cos  v'kV+<^  h(z)  goes  by  Eq.  100  to 

J KdK  e*Kp  cos  + ^ hn(z)  whose  horizontal  asymptote  (p  -»•  °°) 

is  predominantly  along  the  direction  for  which  cos  (cf>  + 6(n))  = 0.) 
Notice  that  as  the  mode  n increases  so  that  £z(n)  decreases,  6(n)-»- 
ir/2  by  Eq.  100,  so  that  $(n)  -*•  0,  by  Eq.  101.  That  is,  the  higher 
modes  contribute  asymptotically  mainly  within  the  downstream  cone 

mm 

<P  < <P(1)  ~ sin  1 , (102) 


which  for  a given  stratified  ocean  varies  inversely  with  V^,  as 
depicted  by  the  stippled  region  in  Fig.  11. 

90° 


60°. 

0 

30° 


0° 

0 1 2 3 4 5 6 

Vb/N2z(n) 

Fig.  1 1 Singular  Wedge  of  Internal-Wave  Wake  of  Horizontal  Azimuth  (0)  as 

a Function  of  Body  Velocity  (Vb)  in  Units  of  Brunt-Vaisala  Frequency 
(N)  Times  Vertical  Scale  of  n*h  Mode  (Cz(n)).  In  general,  ail  singular 
rays  lie  in  the  stippled  region  defined  by  the  lowest  mode  (n  = 1);  for 
the  special  case  of  the  exponentially  stratified  ocean  model,  singular 
rays  for  higher  modes  (n  2,  3,  4 . . .)  are  determined  by  the  dashed 
curves. 
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We  may  illustrate,  and  also  summarize  the  essential  origins 
of  the  singular  behavior,  by  considering  an  exponentially  strati- 
fied ocean  (N  = constant  and  t = D) . In  that  case  we  have  the 
exact  relations, 


Dispersion: 


9 N2k2  N^k2 

“ = ~W  = + kT  * 

z 


Stationarity : to  = kV^  cos  6 ^ , 

Boundary  conditions:  kz  = ~ , (n  = 1,2,...)  , 


which  yield  cos  0Ky  = (Ne/Vjj)//K2e2  + n2^  (cf.  Eq.  39).  Then  for 
p (x  -*•  0) , we  have  for  all  n 


Asymptotically : 


a,  Ne 

cos  6(n)  - 

b 


cf  Eq.  100,  and  by  the  argument  above  the  corresponding 

Singular  directions:  <J>(n)  = sin-1  — ~ , 

V vhu  • 


whence  our  estimate  for  the  downstream  limiting  cone  now  becomes 
exact  (Eq.  102).  In  this  ocean  model,  the  above  results  obtain 
for  all  the  modes  (n  It  1),  as  indicated  by  the  dashed  curves  in 
Fig.  11.  These  preferred  directions  were  earlier  derived  by  R.  I. 
Joseph  from  detailed  calculations  for  the  N = constant,  e = D 
model. 


6.6  ACOUSTIC  WAVES 

The  compressibility  of  seawater  (Table  1)  allows  sound  radi- 
ation to  propagate  with  phase  speed  cs  « 1500  m/s,  so  that 

Vb  <<  cs  (103> 

for  any  normal  ocean  body  motion.  Hence,  in  contrast  to  the  SW 
and  IW  cases  above,  the  stationarity  condition,  Eq.  84,  and  the 
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AW  dispersion  relation,  Eq.  35c,  from  which  w s kV^  and  u ^ kcs, 
respectively,  are  incompatible.  That  is,  in  Fig.  1 the  line  given 
by  Eq.  84  lies  far  below  the  AW  curve,  and  they  cannot  intersect. 
Physically,  no  stationary  wave  system  develops,  because  that  re- 
quires a relative  Doppler  shift  in  the  body  frame  of  unity  (i.e., 
frequency  zero) , whereas  the  maximum  shift  for  AW  is  V^/cs  << 
unity,  by  Eq.  103.  For  the  OM  field  induced  by  AW,  therefore, 
the  body  is  relevant  solely  as  an  approximately  stationary  source 
of  sound,  which  once  emitted  propagates  as  free  AW  through  the 
ocean.  We  thus  neglect  the  body  motion  V^,  entirely  be  Low. 

The  sonomagnetic  pseudowave  emitted  along  with  the  AW  by 
virtue  of  the  interaction  between  acoustic  velocitiec  and  earth's 
magnetic  field  was  discussed  for  plane  waves  in  Subsec.  5.3.  At 
great  distances  from  a radiating  body  such  that  the  AW  is  effec- 
tively plane,  the  results  are  directly  applicable.  In  particular, 
we  may  use  the  simple  OM  estimates  of  Eqs.  56  and  57  for  an  in- 
finite ocf  .1,  the  more  accurate  Eq.  58  and  its  following  discussion 
for  air-sea  and  seabed  boundary  effects,  and  Eqs.  59  for  the  OM 
aerofield  behavior.  Of  course,  the  AW  intensity  level  S = 20  log  6P 
in  Ecs.  57  and  58  is  that  existing  at  the  observation  point.  It  is 
reduced  from  that  at  the  source  position  by  losses,  mainly  the 
transmission  loss  with  proper  allowance  for  channeling  (cf.  the 
discussion  following  Eq.  57).  At  near-to-moderate  distances  from 
the  sound-emitting  body,  however,  the  acoustic  and  sonomagnetic 
waves  cannot  be  considered  simply  plane.  Instead,  more  complex 
fields  are  radiated  as  determined  by  source  and  boundary  geometries. 

In  general,  the  sonomagnetic  field  equations  consist  of  our 
basic  Eqs.  22  for  _h  (with  jg  » jig,  « 0,  and  up  = ,Vij>s)  and  the 

usual  wave  equation  for  the  acoustic  potential  <ps.  For  simple  har- 
monic emissions  of  frequency  w,  these  involve  acoustic  and  magnetic 
propagation  constants, 

k = — and  k = ^-i  you  , (104) 

sc  m 

s 

respectively.  Thus,  the  field  equations  read 

(A  + k 2)h  = yo  curl(H0  x V<J>  ) , (105a) 

div  h = 0 , (105b) 

rs* 

(A  + ks2)$s  = Q , 


(105c) 
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for  any  sonic  source  distribution  Q.  Equations  105  are  subject 
to  the  magnetic  boundary  conditions  (Eqs.  6)  and  appropriate  acous- 
tic ones.  For  given  Q,  the  sonic  field  is  determined  separately 
via  Eq.  105c,  so  we  may  assume  it  to  be  known  and  express  jh  in  terms 
of  4>s  as  follows.  First,  we  derive  a particular  integral  by  writing 
jip  = pa  curKHg  x Vip) , which  identically  satisfies  Eq.  105b  and 
also  satisfies  Eq.~105a  if  (A  + k^H  = $>s.  From  the  similarity 
of  this  to  Eq.  105c,  it  is  easy  to  find  the  solution  \p  = 

(4>s  - " ks2)  where  <J>m  is  the  same  function  as  <{>s  only 

with  ks  ■+  kju  in  it,  i.e.,  (A  + k^2)^  = Q.  Thus,  a particular  in- 
tegral of  Eqs.  105a  and  105b  is 


h 

~P 


-ua  curl[He  * Z(<PS  ~ *n)] 

k 2 
m 


k 2 - 


(106) 


L 


where  ^ For  a monopole  Q (cf.  below),  <j>g  ^ are  just  the 

Green's  functions  for  the  Helmholtz  operators  (A  + k2  m) , respec- 


tively, given  by  Appendix  B,  Eqs.  B.3  and  B.5.  We  can  now  write 
the  complete  solution  for  the  sonomagnetic  field  as 


h + h 


h , 

y 


(within  ocean) 


(above,  below  ocean) 


(107a) 

(107b) 


where  jia  = Vy  with  Ay  = 0 (Subsec.  2.2),  h^  depends  on  seabed  con- 
ductivity, and  hc  is  the  complementary  integral  of  the  homogeneous 
form  of  Eq.  105a,  (A  + k^)!^  = 0,  with  div  ^ = 0 (eq.  105b). 

Thus,  ,hajb,c  can  be  expanded  in  known  functions  - ha  via  y in  har- 
monic functions,  _hc  in  the  transverse  waves  of  the  vector  Helmholtz 
equation,  and  hb  similarly  as  appropriate  (Ref.  47)  - and  finally 
determined  via  the  boundary  conditions. 


For  general  illustration,  we  consider  a monopole  source  in 
a semi-infinite  homogeneous  ocean.  Because  the  superposition  prin- 
ciple is  valid  for  both  acoustics  and  sonomagnetics , the  effects 
of  any  more  general  source  can  be  built  up  as  a sum  of  the  mono- 
pole effects.  Also,  for  a source  whose  submergence  is  shallow  rela- 
tive to  the  ocean  depth,  the  seabed  boundary  is  of  minor  importance 
for  the  OH  field  near  the  body,  so  only  the  air-sea  boundary  need  be 
considered.  Inhomogeneities  of  scale  greater  than  the  source  and/or 


Ref.  47.  P.  M.  Morse  and  H.  Feshbach,  Methods  of  Theoretical 
Physics,  McGraw-Hill,  New  York,  NY,  1953. 


if 


i 

f 

t 
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the  observation  depth  will  also  be  unimportant;  small-scale  inhomo- 
geneity, which  would  greatly  complicate  the  acoustic  and  consequent 
sonomagnetic  fields,  is  for  simplicity  assumed  absent. 

A monopole  emits  AW  with  purely  radially  directed  acoustic 
velocities  up  which  interact  with  the  earth's  field  flg  to  induce 
electrical  current  in  loops  with  a common  axis,  along  JJg.  The  air- 
sea  boundary  may  be  idealized  as  a perfectly  reflecting  plane  for 
the  AW,  equivalent  to  placing  an  image  monopole  above  the  surface 
which  also  induces  current  loops  about  gg*  Thus  physically  one 
sees  that  the  sonomagnetic  field  is  closely  related  to  the  familiar 
electromagnetic  radiation  generated  by  magnetic  dipoles  oriented 
along  .gg,  albeit  in  the  sonomagnetic  case  we  have  to  deal  with  an 
infinite  distribution  of  finite-sized  loops.  The  main  difference 
is  that  the  sonomagnetic  field  will  involve  the  acoustic  propaga- 
tion constant  in  addition  to  the  (low-frequency)  electromagnetic 
one,  ks  and  kjjj  of  Eq.  104.  That  is,  the  sonomagnetic  field  solu- 
tion will  take  a form  analogous  to  the  Sommerfeld  integrals  for 
low-frequency  EM  radiation  from  submerged  radio  antennae,  but  in 
terms  of  both  ks  and  km  instead  of  only  kjj. 

The  sonomagnetic  field  of  the  monopole  is  calculated  in  de- 
tail in  Appendix  B for  the  case  of  vertical  Hg  — appropriate  near 
the  earth's  north  or  south  poles  (indicated  as  + and  -,  respec- 
tively, below).  The  resulting  h field  is  given  via  Eqs.  (107) 

with  hh  irrelevant  here,  hn  determined  by  Eq.  106  from  the  <t  of 

~D  ~P  J n Ts,m 

Appendix  B,  and  aerofield  hg  and  complementary  function  h£  of  the 

form 


h = C (S  ,w)  • eiwt  • I (w,d,p,z)  . (108) 

) L XI  »va  | C 

Here,  C+  is  a constant  depending  on  source  intensity  Sj  (in  db  j[e 

IpPa  at  1 m)  and  frequency  u,  with  |C+I«10  5* 2+^1^20/v2/v2  + vQ2 

in  y - m3  for  typical  ocean  parameteri  (Table  1)  and  v0  w 2 Hz 

(Eq.  55).  .he  I are  Sommerfeld-like  integrals  of  dimensions 
*>■'& ) c 

1/m3  that  ary  with  source  frequency  and  depth  (w,d)  and  the  cylin- 
drical coordinates  (p  horizontal,  z vertical).  These  integrals 
are  not  amenable  to  closed-form  evaluation,  except  for  various 
field-region  approximations  as  given  in  Appendix  B. 

Since  the  outstanding  feature  of  the  sonomagnetic  field 
vis-A-vis  the  sonic  field  is  that  the  former  extends  into  the  air 
above  the  sea  (cf.  the  end  of  Subsec.  5.3),  consider  the  aerofield 
ha  and  in  particular  its  measure  by  a total  field  magnetometer, 
jja  * he  » Jja  * & in  the  polar  regions.  Inserting  parameters  from 
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Table  1 into  Eqs.  B.16  through  B.20  gives  the  following  simple 
approximations.  At  close  ranges  above  a shallow  source  (p,z,d 
« 240/v  and  180/Y\>),  we  have  the  no.npropagating 


a i(i)t 

Near  aerofield:  h • H « a • e 


-10. 6+S . / 2 0 

| a | 10 


ii * + i)A  y 

vR$  y 


> 


(109a) 

(109b) 


where  z(i0)  is  altitude  above  sea  and  R = /p2  + (z  + d)^  is  range 
from  source  (parameters  Si,d,v  * w/2it).  This  is  largest  directly 

overhead  on  the  surface  (p  = z * 0,|a|  ft#  10  10* 6+sl/20/vd) . Still 
nearly  overhead  only  now  at  high  altitudes  (z  » 240/v  and  180/ Yv) , 
|a|  goes  over  from  the  ~1/ (z  + d)2  of  Eq.  109b  to  ~l/zl+  (Eq.  B.17). 

On  the  other  hand,  at  far  horizontal  ranges  but  low  alti- 
tudes (p  >>  240/v  and  180/*^,  z <<  these  values,  and  d < p),  we 
have  the  propagating 


i(o)t-k  p) 

Far  aerofield:  h • H ft#  3 • e S , (110a) 

a**  ci 

a circular  wave  traveling  at  the  speed  u>/ks  = cs  of  sound.  Its 
amplitude  is 


„10-1°-6+Sl/2° 


F6(v) 


in  y 


(110b) 


where  p is  horizontal  range  from  the  source  (Sj,d)  and  Fg(v)  is  de- 
fined in  Eqs.  B.20  and  plotted  in  Fig.  12,  showing  a broad  maximum 
at  v ~ 1 Hz.  The  limiting  frequency  behavior  is  Fs(v)  ft#  v0/v  for 
v >>  Vq  and  2/v/vq  for  v <<  v0,  where  v0  is  the  transition  fre- 
quency («#  2 Hz)  of  Eq.  55.  Therefore,  the  sonomagnetic  pseudowave 
propagating  at  sound  speed  according  to  Eqs.  110  is  greatest  for 
monopoles  ot  infrasonic  frequency  at  deep  submergence  (as  long  as 
d ^ p),  and  of  course  for  high  intensity  Si,  but  falls  off  rapidly 
with  range  p.  For  a monopole  source  Si  = 180  dB  at  v ~ 2 Hz  and 
d ~ 102  m,  Eq.  110b  yields  wave  amplitude  |s|  ~ 2/p2  y. 

More  realistic  acoustic  sources  than  the  monopole  are  in- 
cluded in  the  detailed  sonomagnetic  theory  to  be  presented  in  a 
subsequent  report  of  this  series.  For  example,  the  dipole  sonomag- 
netic pseudowave  can  exceed  Eqs.  110  by  as  much  as  v^/ksd,  so  that 
in  the  numerical  example  above  the  aerofield  amplitude  can  increase 
to  j B j ~ 10/p2  y,  now  independently  of  the  submergence  depth  d. 
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7.  CONCLUSION 


This  survey  substantiates  the  general  conjecture  that  the 
great  majority  of  OM  phenomena  are  reliably  estimated  from  the  rule 
of  thumb: 


h ~ ®M  He  * T U0L  (in  Y)  ’ 


(111) 


for  the  CM  field  magnitude  h in  the  quasi-static  approximation  of 
Eqs.  24,  25,  or  34a.  The  notable  exceptions  are  long  SW  (Subsec. 
5.1.5)  and  infrasonic  AW  (Subsecs.  5.3  and  6.6),  for  which  however 
the  equally  simple  approximation  Eq.  34b  is  applicable. 

The  broad  conclusion  to  be  drawn  from  the  detailed  esti- 
mates of  Secs.  5 and  6 is  that  CM  fields  are  typically  of  the  order 
of  magnitude  of  gammas  or  less.  Despite  this  inherent  weakness  of 
ocean  magnetism,  however,  even  the  smaller  OM  fields  are  consider- 
able relative  to  current  SQUID  (Superconducting  Quantum  Interfer- 
ence Device)  magnetometer  sensitivities  (10-1+  y//Hz  or  better,  cf. 
Refs.  48  and  49  . Thus  in  magnetic  measurements  in  the  ocean  en- 
vironment, all  the  OM  fields  considered  are  potential  contributors. 
Of  course  the  detectability  of  any  one  OM  effect,  in  isolation  from 
the  others  regarded  as  "noise,"  depends  on  the  uniqueness  of  its 
particular  spatial  or  temporal  characteristics. 

The  OM  field  extents  and  their  detailed  spatial-temporal 
configurations  vary  considerably,  with  features  of  special  interest 
from  case  to  case.  The  survey  character  of  this  report  makes  it 
unfeasible  to  summarize  all  of  our  results  here  with  proper  regard 
to  the  assumptions  and  restrictions  made  in  each  case.  It  is  pref- 
erable to  refer  to  the  detailed  discussion  in  the  relevant  subsec- 
tions of  Secs.  5 and  6.  The  glossary  of  major  symbols  given  on  the 
following  pages  will  facilitate  use  of  the  estimation  formulae  we 
have  derived. 


Ref.  48.  J.  Clarke,  "Josephson  junction  detectors,"  Science, 
Vol.  184,  1974,  p.  1235. 

Ref.  49.  J.  Clarke,  "Principles  and  limitations  of  SQUIDs," 
Bull.  Am.  Phys.  Soc.,  Vol.  22-  1977,  p.  86. 
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GLOSSARY  OF  MAJOR  SYMBOLS 

A 

General:  ( ) denotes  vector,  ( ) unit  vector,  ( ) Fourier  amplitude 

A4  AS# 

a wave  amplitude 

(subscripts:  v ■=>  at  frequency  v,  I °=>  internal  wave) 

AW  acoustic  wave 

b subscript  denoting  body  quantity 

c ^ wave  phase  speed 
cg  sound  speed 

d pycnocline  depth,  or  body  depth  (Sec.  6) 

D ocean  depth,  or  body  diameter  (Sec.  6) 

E electric  field 

F (v)  sor.-Mnagnetic  frequency  function  (Fig.  12) 
s 

g gravitational  acceleration 

h induced  magnetic  field 

H total  magnetic  field 

(subscripts:  0 =>  ambient,  e =>  earth,  b *=>  body) 

I magnetic  dip  angle 

IW  internal  wave 

j electric  current  density 

k wave-vector  magnitude 

(subscripts:  s *=>  sound,  m —>  magnetic) 

H acoustic  mode  number 

scale  length  of  IW  mode  n 

L length  scale,  or  body  length  (Sec.  6) 

m magnetic  dipole  moment 

(subscripts:  b **>  body,  hm  ->  hydromagnetic,  0 *=>  combination) 
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buoyancy  frequency  (Brunt-Vaisala) 
ocean  magnetic 


pressure 


radial  distance 


Reynolds  number 

magnetic  Reynolds  number  (Eq.  8) 

radius  of  diffusion  region  for  body  of  velocity  V (Fig.  5) 

acoustic  intensity  level,  or  spin  (Subsec.  4.2),  wave  spectrum 
(Subsec.  5,1.2) 

acoustic  intensity  level  at  1 ir.  from  source  in  db  re  lpPa 

sonomagnetic 

surface  wave 

time  scale  (period  for  waves) 
flow  speed 

hydromagnetie  perturbation  to  the  flow 
body  speed 
Wind  speed 

(subscript:  K =>  in  knots) 
vertical  coordinate 


compressibility 

. 5 

unit  of  magnetism  (=  10  G) 

pycnocline  thickness,  or  dielectric  constant  (Sec.  2) 
shear  vi.cosity 
polar  angle 

angle  between  wave  vector  and  body  velocity 
horizontal  wave  vector  magnitude 
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A wavelength 

y magnetic  permeability 

v circular  frequency 

v0  acoustic  transition  frequency  («*  2 Hz) 

p water  density,  or  horizontal  range 

a electrical  conductivity 

Td  magnetic  diffusion  time 

<p  azimuthal  angle 

<}>s  acoustic  potential 

41  flow  potential 

Y magnetic  scalar  potential 

u angular  frequency 
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Appendix  A 

FORMAL  HYDROMAGNETIC  SOLUTIONS 
IN  QUASI-STATIC  LIMIT 

When  the  hydromagnetic  rate  of  change  (3/3t  ~ 1/T  in  Eqs.  22) 
is  slow  relative  to  the  rate  of  magnetic  diffusion  (1/tj  of  Eq.  8), 
the  describing  Eqs.  22  become 

Ah  « po  curl (H n x u .)  , div  h = 0 , ■ (A. la, A. lb) 


in  the  approximation  of  Eq.  23  and  omitting  a body  field  and/or 
velocity.  Thus  for  slow  temporal  variations,  time  enters  only 
parametrically  in  Eqs.  A.l.  The  first  is  a vector  Poisson  equation 
(Ref.  47,  Chap.  XIV);  using  its  Green's  function  1/ |£  - £ | (i.e., 

A 1/ 1 £ - r | = -4ir6(r  - /))  , one  can  write  at  once  the  particular 
integral 


h (t) 
~P  ~ 


/v  rv 


cur1'  (“o  * JSo^ 


(A.  2) 


By  a parts  integration  (with  uQ  assumed  to  fall  rapidly  enough 
with  distance)  and  from  the  reciprocity  property  (3/3_r  + 3/Sr'  ) 
1/  I.E  “ I = 0»  Eq.  A. 2 becomes 


h 

~P 


(r) 


(A.  3) 


which  evidently  also  satisfies  Eq.  A. lb.  On  the  assumptions  of 
uniform  ambient  field  (Hq  = constc^.t)  and  incompressible  flow 
(div  u0  = 0) , the  use  of  vector  ide-.tities  and  reciprocity  with 
parts  integration  as  above  further  reduces  Eq.  A. 3 to 


h (r)  = (Hn 
~p  ~ 4ir  ~0 


(A.  4) 


corresponuir.6  to  Eq.  26  as  given  in  the  text. 

In  the  quasi-static  approximation  and  on  the  above  assump- 
tions, Eq.  A. 4 thus  gives  the  complete  formal  solution  of  the  hy- 
dromagnetic  Eqs.  22  (w  Eqs.  A.l)  for  an  infinite  medium.  (Recall 
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that  other  hydromagnetic  variables  are  given  in  terms  of  h via 
Eqs.  7.)  In  the  case  of  a bounded  medium,  the  solution  is  modified 
to  satisfy  boundary  conditions  by  adding  to  Eq.  A. 4 a suitable 
linear  combination  of  complementary  solutions  of  the  homogeneous 
form  of  Eqs.  A.l,  i.e.,  the  divergenceless  vector  harmonics  that 
satisfy  the  vector  Laplace  equation  (Ah  • div  h = 0,  cf.  Ref.  47, 
p.  1784). 

In  the  quasi-static  limit,  an  alternative  metnod  of  solution 
is  available  that  determines  h in  two  stages,  instead  of  directly 
via  the  hydromagnetic  Eqs.  5 through  7 (or  22  and  7).  Since  the 
first  stage  of  solution  involves  scalar  Poisson  and  Laplace  equa- 
tions, the  method  can  have  particular  advantage  for  the  applica- 
tion of  boundary  conditions.  The  second  stage  then  determines  h 
via  a Biot-Savart  or  vector  potential  type  integration. 

Thus,  returning  to  Maxwell's  equations  with  3/3 t neglected, 
one  has  from  Eq.  la  that  curl  E = 0,  so  that 

E = -V<)>  (A.  5) 


in  terms  of  an  electrostatic  potential  $.  Taking  the  divergence 
of  Eq.  A. 5 and  using  Eq.  7c  of  the  text  yields 

A<j>  * p divfu  * H)  **  uH.  • curl  un  , (A. 6) 

-w  #s/  U 


the  latter  form  holding  on  present  assumptions  (u«^q»  ^ iJo  = 

constant).  The  approximate  Eq.  A. 6 is  a scalar  Poisson  equation, 
provided  the  unperturbed  flow  is  rotational  (curl  t 0)  and  has 
some  vorticity  along  the  magnetic  firld  (H0  • curl  u0  / 0).  Other- 
wise, notably  for  any  irrotational  flows,  it  becomes  simply  the 
Laplace  equation. 

In  any  event,  with  given  right-hand  side  and  appropriate 
boundary  conditions  in  Eq.  A. 6,  its  solution  determines  E via 
Eq.  A. 5 and  thence  the  current 

j « o(-V$  + uu0  x HJ  (A.  7) 


by  the  modified  Ohm's  law  (Eq.  2a,  and  see  the  discussion  following 
Eq.  4 of  the  text).  Finally,  the  hydromagnetic  field  follows  from 
Eq.  7b  (with  H = constant)  in  the  integral  form 


- r 


h(r) 


r 


(A.  8) 
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One  may  verify  that  Eq.  A. 8 upon  insertion  of  Eq.  A. 7 reduces  to 

the  infinite  medium  Eq.  A. 3,  above,  by  use  of  the  vector  identity 

(with  Eq.  A. 5)  Vd>  x V l/|r  - r I = curl  (E/ 1 r - r*  | ) and  the  inte- 

*W  .'V  *>*  »W  AV  A» 

gral  theorem  P dr'  curl'  (E#/|r-r'|)  = fdo/  x E*  / 1 r - r'  I = 0 for 

field  falling  off  adequately  at  infinity,  and  then  rearranging  the 
remaining  portion  of  Eq.  A. 8 by  use  of  the  reciprocity  property. 
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Appendix  B 

SONOMAGNETISM  INDUCED 
BY  SUBMERGED  ACOUSTIC  MONOPOLE 


| 

| 

$ 


ft 


The  sonomagnetic  field  Eqs.  105  are  solved  by  the  method  of 
Eqs.  106  and  107,  and  the  discussion  following  thereafter,  for  a 
simple  harmonic  acoustic  monopole  submerged  in  a semi-infinite 
homogeneous  ocean,  whose  surface  is  approximated  as  a perfect  sound 
reflector.  Thus,  the  source  term  in  Eq.  105c  is  a monopole  at 
depth  z = -d  and  its  negative  image  in  the  air-sea  surface  (z  = 0) 
at  altitude  z = +d, 

Q = -4irKeia)t  [c(r  + zd)  - 6(r  - zd)  ] . (B.l) 


The  strength  coefficient  K can  be  expressed  via  acoustic  relations 
(Lamb,  Ref.  25)  in  tenns  of  the  source  intensity  level  Sj  (in  dB 
re  1 .iPa)  at  unit  distance  from  the  source  in  the  absence  of  bound- 
ary reflections, 


ic  2 -11+S./20 

K = — — • BP  * 10 
u 


(B.  2) 


where  cg  is  sound  speed,  8 compressibility  of  seawater,  and  P am- 
bient pressure.  Notice  that  the  last  facto*  in  Eq.  B.2  contains 
unit  length,  so  K has  dimensions  of  length3/ time.  The  solution  of 
Eqs.  105c  and  B.l  for  the  acoustic  potential  4>s  is  composed  of 
spherical  waves  from  the  source  and  its  image, 

(io)_-iksr  iwt-iksr/  \ 

" r r ) • (B'3> 

in  which  r is  distance  from  source  and  x from  its  image, 


r = /p2  + (z  + d)2  , r ' - /p2  + (z  - d)2  , (B.4) 


in  cylindrical  coordinates  p,  with  z centered  at  the  a-'r-sea  sur- 
face directly  above  the  source  (at  p = 0,  z = -d) . The  magnetic 
potential  function  used  to  construct  the  particular  integral  h 
in  Eq.  106..  ~ 


P 


m 


= t 


with  k k 


m 


(B.  5) 
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follows  at  once  from  Eq.  B.3.  The  propagation  constants  ks,  km 
are  a s defined  in  Eq.  104.  The  jjp  given  in  Eqs.  B.3  through  B.5 
and  106  is  related  to  the  complementary  and  the  aerofield  .hg 
by  the  boundary  condition  Eq.  10b  at  the  air-sea  interface, 


h = h + h at  z=0, 

»vC 


which  finally  determines  the  complete  solution,  Eqs.  107. 


(B.  6) 


We  consider  specifically  the  case  of  vertical  earth  field. 


5,  - - <«  A 


(B.7) 


approximately  applicable  in  Arctic  and  Antarctic  regions  (+  and  -, 
respectively,  here  and  below).  Our  problem  then  possesses  azi- 
muthal symmetry,  so  that  the  expansions  of  ha  and  hc  indicated  fol- 
lowing Eqs.  107  become 


h = V*  , V = / dq  A(q)  J (qp)eiwt  qZ 

» j*'  u 


in  the  usual  cylindrical  harmonics,  and 


h = curlcurl  z$  , 
~c  ~ 


oo  io)t+/q^-k  ^ z 

4>  = J dq  C(q)  JQ(qp)e 
0 


(B.8a, 
B.  8b) 


(B. 9a) 


(B. 9b) 


in  the  vector-wave  solutions  of  Morse  and  Feshbach  (Ref.  47,  p. 
1823).  Here  A(q),  C(q)  are  coefficients  to  be  determined,  JQ  are 
Bessel  functions,  and  the  positive  root  /q2-^2  is  intended  since 
z S 0 in  the  sea.  As  to  hp,  at  the  air-sea  surface  Eq.  106  with 
Eqs.  B.  3 through  B.5  reduces  for  the  geometry  of  Eq.  B.7  to  the 
form 

UOH2K  /-iksr0  -ikmro\ 

hp  (at  z = 0)  - ±p  3^-2  e1Wt  ^ ^ 1 (B.10) 


5 tr- 

|l  ;5 

r i 

S 6 


where  p is  the  unit  cylindrical  radial  vector  and  r0  = /p2  + d2. 
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Upon  inserting  Eqs.  B.8  through  B.10  into  Eq.  B.6  and  using 
the  Fourier-Bessel  theorem,  one  determines  that 


(B.llb) 

Putting  these  coefficients  into  Eqs.  B.8  and  B.9  gives  the  complete 


solution  to  our  problem  in  terms  of  quadratures  that  are  somewhat 
analogous  in  form  to  the  Sommerfeld  integrals  of  electromagnetic 
theory  for  submerged  radio  antennae.  For  example,  the  aerofield 
scalar  potential  is  from  Eqs.  B.8  and  B.ll: 


Notice  that  the  integral  in  Eq.  B.12a  is  similar  to  that  for  the 
magnetic  aerofield  radiated  from  a submerged,  low-frequency,  verti- 
cal magnetic  dipole  (cf.  Sinha  and  Bhattacharya,  Ref.  50,  Eqs.  3, 

4,  12,  and  13)  but  more  complicated  by  the  appearance  of  us  as  well 
as  um  in  the  integrand.  In  general,  the  integrals  for  the  EM  case 
have  not  been  evaluated  in  closed  form.  Therefore,  a fortiori  in 
the  sonomagnetic  case  we  must  turn  to  approximate  evaluation. 

Consider  that  component  of  the  magnetic  field  in  air  which 
would  be  registered  by  a total  field  magnetometer, 

h = h • H - h • z (B.13) 

ae  ~.jl  ~e  ~ 

for  the  present  case  (Eq.  B.7).  From  Eqs.  B.12a  and  B.8, 


Ref.  50.  A.  K.  Sinha  and  P.  K.  Bhattacharya,  "Vertical  mag- 
netic dipole  buried  inside  a homogeneous  earth,"  Radio  Sci. . Vol. 
1,  1966,  p.  379. 
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-u  d -u  d\ 

S e “ }.  (B. 14) 


m 


).  (B, 


In  terms  of  sonic  and  magnetic  scale  lengths 


o-l  240 

£ = - — w 

s k v 

s 


o - 1 180 


AT 


(B. 15) 


8 


I 


there  are  several  regimes  of  range  coordinates  p,  z,  and  depth  d 
for  which  we  can  give  simple  approximations  to  Eq.  B.14. 


1.  Very  close  above  shallow  source  (p,z,d  << 

^e^61^  * C±  * (km2_ks2)  ' ’ (B*16) 


in  which  C,  • (k  2-k  2),  via  Eqs.  B.12b  ard  B.2,  is 
± m s * 

determined  by  source  strength.  The  regime  of  Eq.  B.16 
includes  the  point  on  the  sea  surface  directly  above 
the  source, 


hae(P=0,z=0)  « eiwt  (?  poHeK)  • , 


provided  now  only  that  the  source  be  shallow  enough 

(d  « £ ,£  ) . 
m s 


2.  High  altitude,  directly  overhead  (z  >>  £m,£s  but  p,  d<^ 
£m»£s) : 


. iwt 

h we  • C.  , 

ae  + \ lk 


-ik  d 

6 \l  e S - e 


m , 


4 


(B. 17) 


which  is  the  same  fall-off  with  height  (l/z1*)  as  the 
"quasi-near"  electromagnetic  field  of  a vertical  mag- 
netic dipole. 


Far  horizontal  ranges  at  low  altitude  (p 


>>  £m»£s» 

jr  • (s.18) 
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This  far-field  circular  wave  reduces  over  measurement  distances 
that  are  sm*  11  compared  to  the  range  (6p  <<  p)  to  a plane  wave 

a where  x = <5p.  The  wave  amplitude  is  from  Eqs.  B.18, 

B.12b,  and  B.2, 

Ih  I w H • BP  • l(f  11+Sl/2°  • F (v)  • -4  y (B. 19) 

1 ae1  e s p2  * 

in  which  the  frequency  behavior  is  contained  in  the  function 

2 (v0/v)  / \-f> 

F (v)  = 1 + x + fl  /I  + x , (B.  20a) 

sv  x \ v v / 


x = A + (v0/v)2  , (B. 20b) 

h^c  2 

v.  = — -x ft#  2 Hz  , (B.20c) 

0 ZTT 

the  last  being  the  transition  frequency  of  Eq.  55,  at  which  k = 

Ul  S 
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